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0\ Abstract 

(N 

We show the existence of the strong inverse correlation between the temporal variations of 

f^ the toroidal component of the magnetic field in the solar tachocline (the bottom of the convective 

r^ zone) and the Earth magnetic field (the Y-component). The possibility that the hypothetical 

_; solar axions, which can transform into photons in external electric or magnetic fields (the inverse 

Q^ Primakoff effect), can be the instrument by which the magnetic field of the Sun convective zone 



modulates the magnetic field of the Earth is considered. 



)—t We propose the axion mechanism of Sun luminosity and "solar dynamo - geodynamo" 

^ connection, where the energy of one of the solar axion fiux components emitted in Ml transition 

i__i in 57pg nuclei is modulated at first by the magnetic field of the solar tachocline zone (due to 

^ the inverse coherent Primakoff effect) and after that is resonantly absorbed in the core of the 

^ Earth, thereby playing the role of the energy modulator of the Earth magnetic field. Within 

t^^ the framework of this mechanism estimations of the strength of the axion coupling to a photon 

1_! {9ay = 7.07 • 10~^^ GeV~^), the axion-nucleon coupling {gan = 3.20 • 10~^), the axion-electron 

■^ coupling [gae = 5.28 • 10""'^"'^) and the axion mass {nia = 17 eV) have been obtained. It is 

.^ also shown that the claimed axion parameters do not contradict any known experimental and 

C^ theoretical model-independent limitations. 

^^ We consider the effect of dark matter in the form of 17 eV axions on the extragalactic back- 

• • ground light. Our treatment is based on theoretical results by Overduin and Wesson (Phys. 

• 5^ Rep. 402 (2004) 267), who described the axion halos as a luminous element of a pressureless 

/\ perfect fluid in the standard Friedman-Robertson- Walker universe basing on the assumption 

S that axions are clustered in Galactic halos with nonzero velocity dispersions. We find that 

the spectral intensity Ix{Xq) of the extragalactic background radiation from decaying axions 

{nia = 17 eV, Ca'yy = 0.02) as a function of the observed wavelength Aq is in good agreement 

with the known experimental data for the near ultraviolet, optical and near infrared bands (1500- 

20000 A). In the framework of such approach it is shown that the present density parameter 

0,a of thermal axions satisfies the inequality 0.12 ^ ila ^ 0.25 and is comparable to the density 

parameter of dark matter. 
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1 Introduction 

In the recent paper by Alessandria et al. the results of CUORE experimental search for ax- 
ions from the solar core from the 14.4 keV Ml ground-state nuclear transition in ^^Fe were 



presented [1]. The detection technique employed a search for a peak in the energy spectrum 
at 14.4 keV when an axion is absorbed by an electron via the axio-electric effect. The cross- 
section for this process is proportional to the photoelectric absorption cross-section for photons. 
In this pilot experiment 43.65 kg ■ d of data were analyzed resulting in a lower bound on the 
Peccei-Quinn energy scale of fa ^ 0.76 • 10^ GeV for the value for the flavor-singlet axial vector 
matrix element of S" = 0.55; bounds are presented in the graph for values 0.15 ^ S" ^ 0.55 
(Fig. 1). With the numbers quoted in the text, the limit on fa translates into the axion mass 
limit rua < 8 eV, significantly more stringent than in the recent results obtained with ^^Fe 
detectors [2, 3] and by the Borexino experiment [4, 5]. 

Despite the fine and elegant experimental implementation of the idea of detecting the solar 
axions through the axio-electric effect in Te02 bolometers (CUORE detection technique [1]), a 
number of fundamental questions regarding the appropriateness of some assumptions used in 
the problem statement arises immediately. 

The first one is rather obvious and lies in the following. Why is the 14.4 keV Ml ground- 
state nuclear transition in solar ^^Fe chosen as the main mechanism of solar axions production 
in CUORE experiment, whereas there are other solar axion production mechanisms discussed 
in scientific literature in detail which also make their respective contributions into the 14.4 keV 
axions flux, such as the so called Primakoff effect (e.g. [2, 3]), bremsstrahlung and the Compton 
process (e.g. [6, 7]) (see Fig. 2)? From the analysis of Fig. 2a, where the spectra of the processes 
under discussion normalized by the corresponding constants are shown, it follows that this 
question is absolutely nontrivial, and the answer depends on the knowledge of the values of 
all these constants simultaneously. In fact, as it will be shown below, the real solar axions 
spectra may look like the ones depicted on Fig. 2b. Therefore the question asked above may be 
reformulated as follows: "What must be the basic physical criterion of the accepted problem 
statement justification, for example, for the experiment on 14.4 keV solar axions detection?" 

In our opinion, one of the most effective ways of establishing such a criterion is the search for 
the models which would describe some experimentally observed phenomena in the framework 
of standard or non-standard solar physics using these properties of axions. If such a model is 
found, then the pivotal estimates of e.g. the axion mass or the upper limits on the axion coupling 
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Figure 1: Expected rate in the axion region as a function of the fa axion constant for different values of the nuclear 
S parameter. The horizontal line indicates the upper limit obtained in CUORE experiment {fa ~ 0.76 • 10^ GeV for 
S = 0.55) [1]. 
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Figure 2: Spectra of solar axions at the ground produced by Primakov effect (.907), Ml ground-state nuclear transition 
in solar ^^Fe {gan)^ bremsstrahlung (gae, green line) and the Compton process (gae, blue line) correspondingly: (a) 
1, gan = 1{S = 0.5), gae = I {S = 0.5); (b) 5„^ = 7.07 • IQ-" GeV, gan = 3.2 • lO'^, g,, = 5.28 • lO"". The 



data by Derbin et al. [7] were used in order to plot the bremsstrahlung and Compton spectra, 
marked with the pink line. 



14.4 keV axions are 



constants to photons {ga-y), nucleons (gan) and electrons (gae), obtained in the framework of 
the given model, may play a role of the main physical justification criterion for the accepted 
problem statement in the 14.4 keV solar axions detection experiment. 

In order to justify such a criterion for the future experiments (e.g. CAST, CUORE, XMASS 
etc.) we decided to create a modified model of the axion mechanism of Sun luminosity^ and 
solar dynamo - geodynamo connection, which had been described in our previous paper [9]. 
The basic idea of such a mechanism, which may be split into two stages for convenience, is 
the following. At the first stage the solar axions fiux variations produced by the previously 
mentioned processes are modulated by the solar tachocline magnetic field variations through 
the inverse Primakoff effect [10]. At the second stage the "modulated" solar axion flux travels 
to the Earth, where its "iron" component containing the 14.4 keV solar axions is resonantly 
absorbed in the iron-nickel core of the Earth. If the energy of the axions supplied to the Earth 
core is enough for generation and maintaining the geomagnetic field, then this process will 
result in a persistent anticorrelation between the variations of the solar magnetic field and the 
geomagnetic field (the Y-component)^ (Fig. 3). This is extremely important, because such effect 
of anticorrelation was discovered recently [9], and has a strong experimental basis. 

It should be added that the solar axion flux modulated by the inverse Primakoff effect in 
the magnetic fleld of the solar tachocline must not only explain the value of solar luminosity, 
but also describe the solar photon spectrum from the Active Sun, which in its turn must be 
equivalent to the data from accumulated observations [14]. 

Thus, the main purpose of the present report was, on the one hand, to develop a modifled 
axion model of the Sun luminosity and solar dynamo - geodynamo connection mechanism; 



^It should be noted here that the axion mechanism of Sun luminosity, which served as a basis for one of the first 
axion mass estimates, was described for the first time in 1978 in the paper [8]. 

^Note that the strong (inverse) correlation between the temporal variations of the magnetic flux in the tachocline 
zone and the Earth magnetic field (the Y-component) are observed only for experimental data obtained at that 
observatories where the temporal variations of declination {6D/6t) or the closely associated east component {6Y/St) are 
directly proportional to the westward drift of magnetic features [11]. This condition is very important for understanding 
of the physical nature of the indicated above correlation since it is known that it is only the motions of the top layers 
of the Earth's core that are responsible for most magnetic variations and, in particular, for the westward drift of 
magnetic features seen on the Earth surface on the decade time scale. Europe and Australia are geographical places, 
where this condition is fulfilled (see Fig. 2 in [11]). For more detailed discussion of this question see below (Section 2). 
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Figure 3: Time evolution of (a) tlie variations of tlie magnetic flux at the bottom (the tachocline zone) of the Sun 
convective zone (see Fig. 7f in Ref. [12]), (b) the geomagnetic field secular variations (the Y-component, nT / year) 
measured at the Eskdalemuir observatory (England) [13]. Curves are smoothed by the sliding intervals in 5 and 11 
years. The pink area is a prediction region. 

and on the other hand, to obtain the consistent estimates for the axion mass and the axion 
coupling constants to photons {ga-f)-, nucleons {gan) and electrons (ffae) through the comparison 
and generalization of the model results and the known experiments including CAST, CUORE 
and XMASS. 



Magnetic field of solar tachocline zone and axion 
mechanism of the solar dynamo — geodynamo 
connection 



It is known that in spite of a long history, the nature of the energy source maintaining a 
convection in the liquid core of the Earth, or more exactly the mechanism of the magnetohy- 
drodynamic dynamo (MHD) generating the magnetic field of the Earth, still has no clear and 
unambiguous physical interpretation [15-19]. The problem is aggravated by the fact that none 
of the candidates for an energy source of the Earth magnetic- field [15] (secular cooling due to 
the heat transfer from the core to the mantle, internal heating by radiogenic isotopes, e.g. ^''K, 
latent heat due to the inner core solidification, compositional buoyancy due to the ejection of 
light elements at the inner core surface) can in principle explain one of the most remarkable 
phenomena in solar-terrestrial physics, which consists in strong (inverse) correlation between 
the temporal variations of the magnetic fiux in the tachocline zone (the bottom of the Sun 
convective zone) [12] and the Earth magnetic field (the Y-component) [13] (Fig. 3). 

At the same time, supposing that the transversal (radial) surface area of tachocline zone, 
through which the magnetic fiux passes, is constant in the first approximation, we can assume 
that magnetic fiux variations also describe the temporal variations of the magnetic field in the 
tachocline zone of the Sun. In this sense, it is obvious that a future candidate for an energy 
source of the Earth magnetic field must not only play the role of a natural trigger of solar- 
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Figure 4: (a) Geomagnetic filed Y-component variations at different observatories [11]; (b) Secular variation of 
declination for 2005-2010 [20] (the closely associated east Y-component of the geomagnetic field), (c) Direct impact 
of the westward drift on the geomagnetic field Y-component in Europe (e.g. Lerwick) and Australia (e.g. Toolangui). 
Here SBy /5t is the magnetic variation seen at magnetic observatories; the dBy/dt term accounts for the effects of 
non-uniform and north-south motions (as well as the effects of magnetic diffusion) , C/^ is a magnitude of the westward 
drift as seen at the Earth surface and dBy /dip is the longitudinal gradient of the magnetic field as seen at the surface. 
(d) Time evolution of the geomagnetic field secular variations (Y-component, nT/year), [13] and the variation of the 
Earth rotation velocity [21] (green line). All curves are smoothed by the sliding intervals in 5 and 11 years. The pink 
area is a prediction region. 

terrestrial connection, but also directly generate the solar-terrestrial magnetic correlation by its 
own participation. 

At this point a question about the physical nature of such correlation arises. Let us turn 
to the concept of the westward drift of the Earth magnetic field. The nondipole part of the 
main field has a characteristic feature - it drifts westward with time. The phenomenon of the 
westward drift was noticed as far back as the XVII century. Each component of the geomagnetic 
filed has its own drift speed with the average westward drift speed 0.2° per year. It means that 
the nondipole field makes one complete revolution around the Earth rotation axis in 1800 years. 
A higher rotation velocity of the mantle in comparison with the outer core is supposed to be the 
physical mechanism of the westward drift. Because of electromagnetic forces, the soHd mantle 
of the Earth is coupled to the core as a whole, and the outer part of the core therefore travels 
westward relative to the mantle, carrying the minor features of the field with it [22]. 

To explain the westward drift of magnetic features we have to distinguish between the drift 
effect and other causes of magnetic variations [11]. The magnetic secular variation can be 
written as: 
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where SBy /Sit is the magnetic variation (the Y-component) seen at magnetic observatories; the 
dBy/dt term accounts for effects of non- uniform and north-south motions (as well as the less 
important effects of magnetic diffusion) , U^p is the magnitude of the westward drift as seen at 
the Earth surface and dBy /dip is the longitudinal gradient of the magnetic field as seen at the 
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Figiire 5: Sketch of Earth magnetic field and westward drift of Earth magnetic field on the core-mantle boundary iu^) 
and at the Earth surface iJJ^). 



surface. 

It is Icnown that most of the early magnetic observatories are located in Europe (Fig. 4a) 
and fortunately, the dBy/d^p term here is large [11] and smoothly varying [11, 23] (Fig. 4b). 
Moreover, the term 5 By /St is dominant and directly proportional to the magnitude of the 
westward drift at the geographical places where the term dBy/dip is large and smoothly varying 
and the term dBy/dt — )■ (for example, Europe and Australia (see Fig. 4c)): 
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Consequently, if the westward drift of the magnetic field on the core-mantle boundary (Fig. 5) 
is caused by the core-mantle coupling, which induces the corresponding westward drift of mag- 
netic feature at the Earth surface (Fig. 5), then 



5By 
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U^ ~ u^- 



(3) 



where u^ is the westward drift of the magnetic field on the core-mantle boundary. 

On the other hand, as far back as 1953 basing on the Bullard's model [22] analysis, Ves- 
tine [24] came to a conclusion that if the core-mantle coupling mechanism exists, it should also 
cause a correlation between the westward drift of the eccentric dipole (the magnetic centre in 
the Earth core) and the variations of the Earth rotation velocity 0,. As the further analysis of 
the magnetic observations and the Earth rotation variations [11] reveals, such kind of correla- 
tion indeed takes place (Fig. 4d), which is an obvious sign of core-mantle coupling mechanism 
existence producing the westward drift of magnetic features at the Earth surface. 

It becomes clear therefore that a question about the physical nature of the strong anticorre- 
lation between the Solar magnetic field and the geomagnetic field (the Y-component) variations 
(Fig. 3) comes to a question: "How can the Sun know about the processes in the Earth liquid 
core?". 

The answer is very simple: it governs them! And it governs the magnetic field in the 
Earth core by means of some unknown interaction carrier! More precisely speaking, the Sun 
governs the processes in the Earth liquid core through some kind of interaction which must 
be transmitted by some unknown particles with their flux controlled (modulated) by the Solar 
magnetic field. 

According to our supposition, these particles may be the axions born primarily inside the 
Sun core and may be converted into 7-quanta in the tachocline magnetic field. This supposition 
is the leading idea of the present paper. 

The fact that the solar-terrestrial magnetic correlation has the undoubtedly fundamental 
importance for evolution of all the geospheres is confirmed by existence of stable and strong 
correlation between temporal variations of the Earth magnetic field, the Earth angular velocity. 



the average ocean level and the number of large earthquakes (with the magnitude M^7), which 
are apparently driven by a common physical cause of unknown nature (see e.g. [9]). 

In this section we consider the hypothetical particles (solar axions) as the main carriers of 
the solar-terrestrial connection, which by virtue of the inverse coherent Primakoff effect can 
transform into photons in external fluctuating electric or magnetic fields [10] . At the same time 
we ground and develop the axion mechanism of solar dynamo - geodynamo connection, where 
the energy of axions, which originate from the Sun core, is modulated at first by the magnetic 
field of the solar tachocline zone (due to the inverse coherent Primakoff effect), and after that 
is resonantly (^^Fe solar axions) absorbed in the iron core of the Earth, thereby playing the 
role of an energy source and a modulator of the Earth magnetic field. Justification of the axion 
mechanism of the Sun luminosity and solar dynamo - geodynamo connection is the goal of the 
current section. 

2.1 Implication from "axion helioscope" technique 
(axion-photon interaction) 

As it is seen from the Earth, the most important astrophysical source of axions is the core 
of the Sun. There, pseudoscalar particles like axions would be continuously produced in the 
fluctuating electric and magnetic fields of the plasma via their coupling to two photons (the 
Primakoff effect [10]). After production the axions would freely stream out of the Sun without 
any further interaction. The resulting differential solar axion flux on the Earth would be [25, 26] 
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where E is in keV and giQ = ga^y/iW'^^ GeV~^). 

The spectral energy of the axions (4) follows the thermal energy distribution between 1 
and 100 keV, which peaks at ss3 keV and the average energy (£'„) = 4.2 keV. To be able to 
compare the expected axion flux in a specific energy range with available data, by integrating 
the spectrum (4) over the energy range of 1 to 100 keV we find the solar axion flux at the Earth 
to be 
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In the case of the coherent Primakoff effect the number of photons leaving the magnetic field 



towards the detector is determined by the probability Pa- 
" observable" photon inside the magnetic field [27] 
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where B is the strength of the transverse magnetic field along the axion path, L is the path 
length traveled by the axion in the magnetic region, / = 27r/q is the oscillation length, F = A~^ 
is the absorption coefficient for the X-rays in the medium, A is the absorption length for the 
X-rays in the medium and the longitudinal momentum difference q between the axion and the 
X-rays energy E^ = Ea is 



with the effective photon mass 
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where m-y = rua is the axion mass, a is the fine-structure constant, Ue is the number of electrons 
in the medium, rrie is the electron mass, Z is the atomic number of the buffer medium, A is 
atomic mass of the medium and its density p in g/crrfi . 

On the other hand, it is known that the axion is a neutral pseudoscalar particle that was in- 
troduced in the particle theory to explain the absence of CP violation in strong interactions [28- 
30]. The most natural solution to the CP- violation problem was obtained by introducing a new 
chiral symmetry, known as Peccei-Quinn (PQ) symmetry [1], the spontaneous breakdown of 
which at the energy fa fully compensates the CP-nonivariant term in the QCD Lagrangian and 
leads to the appearance of the axion [29, 30]. The axion is not massless because the chiral U(l) 
PQ-symmetry is anomalous. As a result, the axion gets a mass of the order [31] 

QCD ,(.. 

ma ~ —, — , (9) 

Ja 

where Aqcd is the confining QCD-scale and fa is the energy scale associated with the breakdown 
of the U(l) PQ symmetry. 

At the same time it is necessary to mention the axion mass estimates obtained in the 
framework of the so-called invisible axion models (KSVZ [32, 33] and DFSZ [34, 35]), which 
restrict the allowed range for fa, or equivalently the range for the axion mass 

10^ GeV 
ma = 6-^^^ eV, (10) 

Ja 

2.2 Axion conversion in the Sun magnetic field and the 
plasma mass of photon 

Let us consider the modulation of the axion flux emerging from the Sun core and passing through 
the solar tachocline region (ST) located at the base of the solar convective zone (Fig. 6c, d). As 
is known [14, 36], the equatorial thickness of ST, where the toroidal magnetic field B ~ 10-^50 T 
dominates [36-38], attains Lst ~ 0.039Rs (where Rs = 6.96 • 10^ m [39] is the Sun radius). At 
the same time the values of pressure, temperature and density for the ST are Pst ~ 6.0-10^^ Pa, 
TsT ~ 2.0 • 10^ K and p ~ 0.2 g • cm~^, respectively. 

To estimate the plasma mass of a photon m-y in the hydrogen-helium medium of ST it is 
possible, without loss of generality, to use the modified Eq. (8) in the form [25] 



m,{eV) = ma = Jo.Ol^f^'^ - 25 eV, (11) 

where we use the corresponding parameters Pst ~ 6.0 ■ 10^^ Pa and Tst ~ 2.0 • 10^ K for the 
hydrogen-helium medium of ST obtained by Bahcall &: Pinsonneault for the standard model of 
the Sun [41]. 

Thus, the axion mass in the standard model of the Sun is ~25 eV. However, it will be shown 
below that in the framework of the axion mechanism of Sun luminosity the axion mass will be 
different, since the total energy balance of the Sun is not violated, but indicates a substantial 
change in radiation transport through the radiative zone and the convective zone with respect to 
the standard model of the Sun. The energy portion of the axion-independent radiation transport 
is rather small here and equals to ~ 0.015Asun (see (19)). Since the total energy balance of 
the Sun is not violated in the axion model, one may suppose that the basic parameters of the 
solar core - the region of the energy generation - remain approximately the same in both the 
standard and the "axion" models. Meanwhile, the thermodynamic parameters (temperature, 
pressure, plasma density, electron density etc.) outside the solar core (between the core and 
photosphere) are substantially smaller in the axion model as compared to the corresponding 
parameters in the standard model of the Sun. 

It should be noted here that the calculation of these parameters for the axion model is 
a rather nontrivial task, which, because of its complexity, will be performed in a separate 

9 




Figure 6: Examples of simulation of the periodic alternation of layers (the zonal flow) in the convective structures of 
the Earth outer core (a, b [40]) and the convective zone of the Sun (c, d). a) View from the north. Isosurfaces of the 
axial vorticity, W2, are shown in red (cj^ — 0.4) and green (wz = —0.4) to illustrate the sheet plumes. Each line forms 
a closed ring, indicating that the flow is nearly purely westward; b) Same as a), but viewed from a different angle; c) 
The section AA' along one of the alternate convective layers of the Sun. In the tachocline axions are converted into 
7-quanta (see the inset in d)) channeling in the green area. In the photosphere 7-quanta are scattering due to the 
Compton effect, d) Same as c), but viewed from a different angle (see the section A A' in c)). Alternate layers in the 
convective zone, where layers in which the channeling takes place are shown in green, are also presented. Blue points 
(in the upper convective zone) and crosses (in the tachocline) show the direction of output and input of the magnetic 
field in the convective zone of the Sun. 

publication. For this reason from now on let us use the "experimental" value of the axion mass 
found during the extragalactic background light investigation (see Section 5). 



rua ~ 17 eV. 
Taking into account (11) and (10) it is easy to derive the value of the energy /„: 



(12) 



/a = 0.353 -10^ GeV 



(13) 



Now we make an important assumption that the axion mass is equal to the plasma mass 
of a photon, i.e., m^ = rUa ~ 17 eV. It is obvious, that by virtue of Eq. (7) g — t- 0, whence it 
follows that the oscillation length / becomes an infinite quantity, i.e. / = 2TT/q — t- c«. However, 
taking into account that in this case the absorption length A is about 0.1 m [14], we have 
TLsT -^ 00. This means that according to Eq. (6), the intensity of expected conversion of 
axions into 7-quanta is practically equal to zero in this case. 

At the same time, there is a reason to believe (see [14] and Refs. therein) that the conversion 
of axions into 7-quanta indeed takes place, and strangely enough, this process goes on quite 
effectively. For example, the reconstructed solar photon spectrum below 10 keV from the Active 
Sun (Fig. 7b) is well described by the sum of secondary Compton spectra obtained e.g. by the 
simulation of 7-quanta passage (regenerated from the solar axion spectrum in the tachocline 
zone of the Sun (Fig. 7b)) through the areas of the solar photosphere of different thickness but 
equal density, layers with the thickness of 64 g/cm? and 16 g/cm?. 

In other words, despite the fact that the coherent axion-photon conversion by the Primakoff 
effect is impossible due to the small absorption length for 7-quanta (/ ^> 1) in the medium (see 
Eq. (6) r = 1/A — )■ 00), there is a good agreement between the relative theoretical 7-quantum 
spectra generated by solar axions and experimental photon energy spectra detected close to the 
Sun surface in the period of its active phase (see Fig. 7). The additional account taken of the 
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Figure 7: Reconstructed solar photon spectrum below 10 keV from the Active (flaring) Sun (the black line) from 
accumulated observations [42] (adapted from [14]). The dashed line is the converted solar axion spectrum. Two 
degraded spectra due to multiple Compton scattering are also shown for column densities above the initial conversion 
place of 64 g/cm? , 16 g/crn?. The pink dotted line represents the initial Primakoff axion spectrum. Note that the 
Geant4 code photon threshold is at 1 keV and therefore the turndown around ^-1 keV is an artifact. 

bremsstrahlung in the photosphere will surely enhance the quality of the theoretical description 
of the experimenal solar photon spectrum substantially. 

At the same time it is necessary to note that attempts to match the absolute values of these 
spectra did not succeed so far [14]. It was mainly associated with the absence of a wish to make 
efforts, since it was absolutely unknown how can the 7-quanta spectra generated by solar axions 
in the tachocline be transported in a "virgin" , i.e. unchanged, form through the convective zone 
up to the Solar photosphere. 

To overcome the problem of the small absorption length for 7-quanta and to reach a reso- 
nance in Eq. (6) it is necessary for the refractive gas, in which the axion-photon oscillation is 
studied, to have a zero refractive index [43]. It appears that in order to satisfy this condition 
it is not necessary to use the so-called metamaterials [44] with the negative permittivity [e) 
and magnetic permeability (^) or results of the Pendry super lens theory [45], which are not 
practically realized in nature'^. Taking into account the known difficulties [51] induced by the 
so-called problem of electromagnetically induced transparency for X-rays and the recent signif- 
icant advances in this field [52, 53], let us consider two (possibly related) alternative ways of 
solving the problem of "unperturbed" 7-quanta spectra transfer through the Solar convective 
zone. 

2.3 Channeling of 7-quanta in periodical structure 

We can use the results from papers [54, 55], where the possibility of the electromagnetic X- 
radiation in a microwave range channeling in a multi-layered metal-dielectric structure is theo- 
retically and experimentally shown. 

^Though it should be noted that the metamaterial technology is frequently used nowadays for laboratory simulations 
of some celestial mechanics and cosmology phenomena [46-49]. And the "..."artificial atoms" used as building blocks 
in metamaterial design offer much more freedom in constructing analogues of various exotic spacetime metrics, such 
as black holes, wormholes, spinning cosmic strings, and even the metric of Big Bang itself. Explosive development 
of this field promises new insights into the fabric of spacetime, which cannot be gleaned from any other terrestrial 
experiments" ([50] and Refs. therein). 
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As it is stated in Appendix A in detail, the essence of the electromagnetic X-ray channeling 
in long-period media lies in a fact that the rays are reflected from the layers of higher electron 
density when propagating at small angles to these layers. It leads to a non-uniform intensity 
distribution over the cross-sectional plane because of the rays concentration within the " chan- 
nels" - the layers with lower electron density. It decreases the absorption substantially and 
makes it possible for the rays to penetrate much deeper into the sample along the layers than 
in the case of an arbitrary angle of arrival. 

According to [54], the intensity J(x) of the photons (see Fig. A. 2 and (A. 26)) passed through 
a sample of a thickness x, may be written in the form 



J{x) = Joexp(-o-x) = Joexp ( x) ■Q{a,yo,x), (14) 

cos a 



where 



Q{a,yo,x] 

■ exp 



x^/32^ (l _ ^i^rll' 



^xi^/32^ {l _ EM.' 



at q>l, 

(15) 
at q < I. 



with the same notation used in expressions (A.26)-(A.27). 

Although we give a complete analysis of the Eqs. (14) and (15) in A, let us make a short 
remark regarding the physical nature of these equations. Here the multiplier Jq exp(— xx/ cos q) 
in (14) corresponds to the case of 7-quanta propagation through a homogeneous medium with 
the electron density A'^g and the absorption coefficient xo- The additional multiplier Q{a,yo,x) 
characterizes the influence of the medium layering. 

As is shown in A, the condition Q(a,yo,x) — >• 1 is theoretically feasible for a majority of 
the multilayer metal-dielectric structures [54, 56], which are an effective emulator of a plasma 
medium (Fig. 6). This condition is obviously necessary, but not sufficient. The layers with 
ultralow, if not with "quasi-zero" density, are also required for the ideal photon channeling. 
Such layers suppress the photon absorption processes almost completely, i.e. minimize the 
effect of the multiplier Joexp(— xoa^/coso) in (14). 

Surprisingly enough, it turns out that such long-period (in terms of density) media with one 
of the two alternating media having almost zero density can take place, and not only in plasmas 
in general, but straight in the convective zone of the Sun. Here we generally mean the so-called 
magnetic flux tubes, the properties of which are examined below (see Appendix B for details). 

2.4 Channeling of 7-quanta along the magnetic flux tubes 
(waveguides) in Solar convective zone 

The idea of the energy flow channeled along a fanning magnetic field has been suggested for the 
first time by Hoyle [57] as an explanation for darkness of umbra of sunspots. It was incorporated 
in a simple sunspot model by Chitre [58]. Zwaan [59] extended this suggestion to smaller flux 
tubes to explain the dark pores and the bright faculae as well. Summarizing the research of the 
convective zone magnetic fields in the form of the isolated flux tubes, Spruit and Roberts [60] 
suggested a simple mathematical model for the behavior of thin magnetic flux tubes, dealing 
with the nature of the solar cycle, the sunspot structure, the origin of spicules and the source of 
mechanical heating in the solar atmosphere. In this model, the so-called thin tube approximation 
is used (see [60] and Refs. therein), i.e. the field is conceived to exist in the form of slender 
bundles of field lines (flux tubes) embedded in a field- free fluid. Mechanical equilibrium between 
the tube and its surrounding is ensured by the reduction of the gas pressure inside the tube, 
which compensates the force exerted by the magnetic field. In our opinion, this is exactly the 
kind of mechanism Parker [61] was thinking about when he wrote about the problem of fiux 
emergence: "Once the field has been amplified by the dynamo, it needs to be released into the 
convection zone by some mechanism, where it can be transported to the surface by magnetic 
buoyancy" [62]. 
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In order to understand magnetic buoyancy, let us consider an isolated horizontal flux tube 
in pressure equilibrium with its non-magnetic surroundings, so that 

Pint + 7, = Pext, (16) 

2/Uo 

where pint and pext are the internal and external gas pressures respectively and ^uq is the magnetic 
permeability of the medium, B denotes the uniform field strength in the flux tube. If the 
internal and external temperatures are equal so that Tint = Text (thermal equilibrium), then 
since pint < Pext, the gas in the tube is less dense than its surrounding (pint < Pext), implying 
that the tube will rise under the influence of gravity. 

In spite of the obvious, though turned out to be surmountable, difficulties of expression (18) 
application to the real problems, it was shown (see [60] and Refs. therein) that strong buoyancy 
forces act in magnetic flux tubes of the required field strength (10^-10^ G [63]). Under their 
influence tubes either float to the surface as a whole (e.g. Fig.l in [64]) or they form loops of 
which the tops break through the surface (e.g. Fig.l in [59]) and lower parts descend to the 
bottom of the convective zone, i.e. to the overshoot tachocline zone. The convective zone, being 
unstable, enhances this process [65, 66]. Small tubes take longer to erupt through the surface 
because they feel stronger drag forces. It is interesting to note here that the phenomenon of 
the drag force which raises the magnetic flux tubes to the convective surface with the speeds 
about 0.3-0.6 km/s was discovered in direct experiments using the method of time-distance 
helioseismology [67]. Detailed calculations of the process [68] show that even a tube with the 
size of a very small spot, if located within the convective zone, will erupt in less than two years. 
Yet, according to [68], the horizontal fields are needed in the overshoot tachocline zone, which 
survive for about 11 yr, in order to produce an activity cycle. 

A simplified scenario of magnetic fiux tubes (MET) birth and space-time evolution (Fig. 8a) 
may be presented as follows. MFT is born in the overshoot tachocline zone (Fig. 8d) and rises 
up to the convective zone surface without separation from the tachocline (the anchoring effect), 
where it forms the sunspot (Fig. 8b) or other kinds of active solar regions when intersecting the 
photosphere. There are more fine details of MFT physics expounded in overviews by Hassan [62] 
and Fisher [64], where certain fundamental questions, which need to be addressed to understand 
the basic nature of magnetic activity, are discussed in detail: How is the magnetic field generated, 
maintained and dispersed? What are its properties such as structure, strength, geometry? What 
are the dynamical processes associated with magnetic fields? What role do magnetic fields 
play in energy transport? 

Dwelling on the last extremely important question associated with the energy transport, 
let us note that it is known that the thin magnetic flux tubes can support longitudinal (also 
called sausage), transverse (also called kink), torsional (also called torsional Alfven), and fluting 
modes (e.g. [72-76]); for the tube modes supported by wide magnetic flux tubes see Roberts and 
Ulmschneider [75]. Focusing on the longitudinal tube waves known to be an important heating 
agent of solar magnetic regions, it is necessary to mention the recent papers by Fawzy [77], 
which showed that the longitudinal flux tube waves are identifled as insufficient to heat the 
solar transition region and corona in agreement with previous studies [78]. 

In other words, the problem of generation (the source) and transport of energy by magnetic 
flux tubes remains unsolved in spite of its key role in physics of various types of solar active 
regions. 

Interestingly, this problem may be solved in the natural way in the framework of the " ax- 
ion" model of the Sun. As it is shown in Appendix B, the inner pressure, temperature and 
matter density decrease rapidly in a magnetic tube "growing" between the tachocline and the 
photosphere. The analysis of these parameters evolution within the equation of the growing 
magnetic flux tube medium state not only gives the ultralow values for them, but also leads to 
the so-called hydrostatic condition of an ideal (without absorption) 7-quanta channeling inside 
the thin magnetic flux tubes 



13 







Figure 8: (a) Vertical cut through an active region iUustrating the connection between a sunspot at the surface and 
its origins in the toroidal field layer at the base of the convection zone. Horizontal fields are stored at the base of 
the convection zone (the overshoot tachocline zone) during the cycle. Active regions form from sections brought up 
by buoyancy (one is shown in the process of rising). After the eruption through the solar surface a nearly potential 
field is set up in the atmosphere (broken lines), connecting to the base of the convective zone via almost vertical flux 
tube. Hypothetical small scale structure of a sunspot is shown in the inset (adopted from Spruit [69] and Spruit and 
Roberts [60]). 

(b) Detection of emerging sunspot regions in the solar interior [67]. Acoustic ray paths with lower turning points 
between 42 and 75 Mm (1 Mm=1000 km) are crossing the region of the emerging flux. For simplicity, only four out 
of a total of 31 ray paths used in this study (the time-distance helioseismology experiment) are shown here. Adopted 
from [67]. 

(c) Emerging and anchoring of stable flux tubes in the overshoot tachocline zone, and its time-evolution in the 
convective zone. Adopted from [70]. 

(d) Vector magnetogram of the white light image of a sunspot (taken with SOT on a board of the Hinode satellite - 
see inset) showing the direction of the magnetic field and its strength (the length of the bar) in red. The movie shows 
the evolution in the photospheric fields that has led to an X class fiare in the lower part of the active region. Adopted 
from [71]. 



Pext 



\Bl_ 
2/xo 



(17) 



which is well satisfied for the "axion" model of the Sun, according to estimations in Appendix B. 
It means that such thin magnetic flux tubes are the ideal 7-quanta waveguides, which reveal the 
essence of the unique energy transport mechanism between the tachocline and the photosphere. 
As a matter of fact, the phenomenon of 7-quanta channeling along the magnetic flux tubes 
not only makes it possible to solve a problem of the energy transport to the photosphere, but may 
also be a basis for solving other important and critical problems in solar physics. If we assume 
that the vertically oriented thin magnetic flux tubes play the role of waveguides for 7-quanta 
produced in the tachocline via the axion mechanism of Sun luminosity, virtually all known 
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Figure 9: (a) Growth rates for magnetic shear instabihties are plotted as a function of the initial latitude (vertical 
axes) and the field strength (horizontal axes) of a toroidal band. Shaded areas indicate instability in 0.1-lOOT band 
(gray) and 200-400T band (green). Contour lines represent m = 1 and m — 2 symmetric (S) and antisymmetric (A) 
modes as indicated. The non-dimensional model is normalized in such a way that the growth rate of 0.01 corresponds 
to an e-folding growth time of 1 year. The parameter s is the fractional angular velocity contrast between equator 
and pole and the reduced gravity G (adopted from [79]). In addition, a hidden part of the "latitude - magnetic field 
in overshoot tachocline zone" dependence, which was missing on the original plot (Fig. 11 in [79]), is plotted to the 
right of the dashed line. 

(b) Solar images at photon energies from 250 eV up to a few keV from the Japanese X-ray telescope Yohkoh (1991- 
2001) (adapted from [14]). The following shows solar X-ray activity during the last maximum of the 11-year solar 
cycle. 

anomalies of experimental data interpretation in physics of active solar regions, helioseismology 
and solar neutrino are withdrawn. Since this assumption needs to be substantiated, let us 
describe our phenomenology, consequences and experimental proofs of this hypothesis below in 
short. 

First of all, if one takes into account the sufficiently strong magnetic field in the balance 
equation of (16) type, it becomes clear that the vertically oriented thin magnetic flux tubes may 
serve as the X-ray waveguides for the radiation originating from the overshoot tachochne zone 
because of the high magnetic pressure. Naturally, in this case the X-ray spectrum coincides with 
the observed Solar X-ray spectrum. All these facts, i.e. the strong magnetic field ~200-400T 
(Fig. 9), X-rays (Fig. 8b) and their spectrum (e.g. Fig. 7 in [71]) in active solar regions, are in 
good agreement with observational data. 

Second , it clears up a way to the solution of the known problem associated with the over- 
shoot tachocline anomaly (Fig. 10) which arises when interpreting the helioseismology and solar 
abundances data. And here is why. 

It is known [80] that the problem comes from the attempts to improve agreement between 
solar models with low heavy-element abundances and seismic inference. The low-metallicity 
models that have the least disagreement with seismic data require changing all input physics 
to stellar models beyond their acceptable ranges. Let us consider the way it happens in the 
framework of a solar model built upon the axion mechanism of Sun luminosity. 

Since solar luminosity is determined by the 7-quanta born in the tachocline in the framework 
of the axion mechanism, it is clear that an old heat flux transport mechanism (from the radiative 
zone to the overshoot) by radiation, used in the standard model of the Sun, should be highly 
depressed because the major part of the radiation is converted into axions in the core of the 
Sun [25] and does not get into the radiative interior. It is easy to see from the traditional state- 
ment of the problem involving helioseismology and solar abundances as described by Basu [80] 
that this is one of the main and fundamental differences from the standard model of the Sun: 
"The most easily detectable effect of the reduction of heavy-element abundances is a change 
in the position of the base of the convection zone. The temperature gradient in the radiative 
interior is determined by opacity, and hence, its structure is affected by the heavy-element abun- 
dances. The base of the convection zone occurs at a point where opacity is just small enough to 
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Figure 10: The relative sound-speed (the panel a) and density differences (the panel b) between the Sun and the model 
constructed with the AGS05 abundances [81]. For comparison we also show the results for the model constructed with 
the GS98 abundances. MDl 360 day data have been used for the inversions. The overshoot tachocline anomaly is 
highlighted with green. Note: GS98 is a solar model with the solar heavy-element mixture Z/X = 0.0245; AGS05 is 
a solar model with low heavy-element mixture Z/X ~ 0.0122. Adopted from [80]. 

allow the entire heat flux to he transported by radiation, and thus the location of this point 
depends on the abundance of those heavy elements that are the predominant sources of 
opacity in that region. If these abundances are reduced, opacity reduces, and the depth of the 
convection zone also reduces. Since the depth of the convection zone has been measured very 
accurately, it is the most sensitive indicator of opacity or heavy-element abundances" . 

The axion mechanism of Sun luminosity implies that because of a virtually complete trans- 
parency of the magnetic flux tubes for 7-radiation there are no reasons for moving the location 
of the center of the tachocline "by hand", since the radiative opacity determined by the effect 
of heavy-element abundance loses its impact on the location of this point and, consequently, its 
significance, because of almost total suppression of the radiative heat flux transport mechanism 
itself in this case. In other words, the effect of absolute magnetic flux tubes transparency for 
the 7-radiation is dominant in the overshoot tachocline zone and levels the influence of radiative 
opacity. As a consequence, the value of the temperature gradient in the radiative interior is al- 
most entirely free from the strict limit introduced by opacity which is still affected by the heavy 
element abundances, but not as dramatically as it is in several other solar models - standard 
and nonstandard - that have been published recently [80]. The latter opens up a possibility to 
build a new standard solar model on the basis of the axion mechanism of Sun luminosity which 
may become a key to the solar abundance problem solution. 

Third , let us consider the axion mechanism of Sun luminosity compatibility with the nuclear 
energy generation pathways in the solar core and solar neutrino fluxes generation. 

The axion mechanism of Sun luminosity is compatible with the standard nuclear energy 
generation pathways scheme and does not disturb the known values for solar neutrino fluxes, 
since the " invisible" axion losses almost do not change the Sun energy balance in our model (see 
Section 2.5 below), and therefore do not introduce any problems related to energy- producing 
regions (i.e. the solar core). 

It actually means that introducing the axion mechanism of Sun luminosity in the framework 
of the standard model of the Sun leads to such value of axion losses which does not contradict the 
Gondolo-Raffelt limit on the "invisible" axion and Sun luminosities ratio, A^J^'"^^ / Asun ^0.1 [82], 
for which a good coincidence between the theoretical values and experimental data of modern 
helioseismological and solar neutrino experiments is still observed [82, 83]. 
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And forth , if the vertically oriented thin magnetic flux tubes in the convective zone play 
a role of the waveguides for the 7-quanta born in the tachocline with total luminosity equal 
to that of the Sun, what are the nature and the power of the energy source maintaining the 
convective processes on the Sun? 

In order to find it out, let us assume that this source is the radiative zone and perform 
an estimation of its power basing on the magnetic field Bqt in the overshoot tachocline zone 
dependence on the total ohmic dissipation D'^^^'^ in the convective zone [84]. 



D"CZ = y - r ^ ^Ot) dV OC j^Ernag, (18) 

where Emag is the magnetic energy of the field which could be possibly maintained by the 
currents that produce the ohmic dissipation of the solar dynamo [85] . 

It is easy to show [84] that the expression (18) may be written down in the following form: 

j^oh^ic ^ iL^Bl « 0.015As.„, (19) 

where D^^^'^ = 0.015A5„„ is the heat power of the radiative zone near the border of the 
overshoot tachocline zone, equal to the uncertainty of the known Solar luminosity [76, 86, 87]; 
r] ~ 10^ crv? I s is the magnetic diffusivity [88, 89], Vcz is the volume of the Sun convective 
zone; /x ~ 1 is the permeability; Bqt = 400 T is the magnetic field in the overshoot tachocline 
zone; Hp = 6.5 • 10^ km is the pressure scale height ■* [59, 87] 

The approximate equality implies that not only the total solar energy balance is preserved 
in the framework of the axion mechanism of Sun luminosity, but also that the temperature 
transport changes substantially with respect to the standard model of the Sun. This is because of 
the fact that the old heat flux transport mechanism (from the radiative zone to the overshoot) by 
radiation, used in the standard model of the Sun, is highly depressed because the majority of the 
radiation is converted into axions in the core of the Sun and therefore does not reach the radiative 
interior. At the same time, this change may not seem so dramatic, since almost all known 
anomalies of the experimental data interpretation on the active solar regions, helioseismology 
and solar neutrino may be leveled as it was noted above. 

And, finally, turning back to the possible mechanisms of 7-quanta channeling in a periodical 
structure and along the magnetic flux tubes in the Sun convective zone, one may suppose with 
confidence that they are not only physically compatible, but may turn out to be just two different 
versions of the same mechanism, which naturally manifests itself, for example, in the so-called 
hexagonal magnetoconvection kinetics (see e.g. [93]). 

This means, in its turn, that the absorption length A for photons in such a medium (see (6)) 
will become considerably greater than the thickness of the overshoot tachocline zone, i.e., A S> 
Lqtz- At the same time, it is obvious that F = A~^ — )• 0, whence a necessary condition 
TLqtz — ^ follows. In the particular (non-coherent) case in which the magnetic field where 
axions are converted into photons is under vacuum (T — t- 0, m^ — t- 0), equation (6) becomes 

F„,= (?=^^^f^)'sm^(?^)/(^)' (20) 

where q = rn^/2Ea (see (7)). 

Obviously, in coherent case (7 — )• 0, regardless of the 7-quanta channeling mechanism type, 
the probability (20) for an axion to be converted back to an "observable" photon inside the 
magnetic field may be expressed in the following simple form 



^A larger value of the pressure scale height is a consequence of the fact that the rigidity [90] of the interior can be 
provided only by the large-scale magnetic field (cf. Mestel & Weiss [91]; Gough & Mclntyre [92]) that the tachocline 
provides the interface in which radial field lines might connect the convection zone with the radiative interior only 
near the latitudes at which there is essentially no radial shear. 
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P..-(^^^^^%^)', (21) 

where Bqt is the mean value of the magnetic field in the overshoot tachochne zone with the 
effective thickness Lqt- A value for Lqt was chosen basing on the analysis of the following 
data set. The most well known results obtained by Charbonneau et al. [94] yield a tachocline 
thickness of At/Rs = 0.039 ±0.013 at the equator and At/Rs = 0.042 ±0.013 at the latitude of 
60°, suggesting that the tachocline may get somewhat wider at high latitudes but that the result 
is not statistically significant. On the other hand, Basu and Antia [95] argue for the statistically 
significant increase in the tachocline thickness with the latitude, from At/Rs ~ 0.016 at the 
equator to At/Rs ~ 0.038 at latitudes of 60° (when the width is defined as in [95]). Furthermore, 
they suggest that the variation may not be smooth; there may be a sharp transition from a 
narrow tachocline at low latitudes to a wider tachocline at high latitudes, possibly associated 
with the sign of the radial angular velocity gradient which reverses at the latitude of ~ 35°. 
Other estimates for the width of the tachocline range from O.Oli?^ to 0.09i?5 (Kosovichev [96], 
Basu [97], Corbard et al. [98], Elhott and Gough [99], Basu and Antia [100]). 

Taking into account that, first, the tachocline is a transition layer between two distinct 
rotational regimes (the differrentially rotating solar envelope and the radiative interior) where 
the rotation is uniform, second, the maximum estimate of the tachocline thickness reaches 
0.09i?5 and third, the thickness of the overshoot tachocline zone is somewhat larger than that 
of the tachocline, we took the value of Lqt equal to O.I-R5. 

Then using Eq. (21) and the parameters of the magnetic field, it is possible to write down 
the expression for the solar axion flux"' probability at the Earth as 

P^ =lf 9jn \UBotV( Lot \\ , . 

"^^ 4 V7.07-10-" GeF-i/ V400 Ty* \7.2b-Wm) ^ ' 

where the value of the magnetic field Bqt = 400 T (cf . [36-38] ) was chosen so that it satisfied 
the "experimental" estimates of 200-400 T (see. Fig. 9) and induced via (22) such a value of the 
axion-photon coupling constant {ga-y = 7.07 • 10~^^ GeV~^) that would in its turn be strictly 
consistent with the known and very important limits (86)-(87) taking into account (13). More 
detailed justification of such self-consistent choice will be given in Section 5. 

It is necessary to make a deviation concerning some important features of the oscillation 
length (/ = 27r/g) here. It is known that in order to maintain the maximum conversion proba- 
bility, i.e. zero momentum transfer (^ — )• 0), the axion and photon fields, put into some medium 
(m-y = ma), need to remain in phase over the length of the magnetic field. This coherence 
condition is met when qL ^ vr, and along with (7) lets one obtain the following remarkable 
relation [14] between the medium density variations and axion mass variations for the coherent 
case of g — )• 0, i.e. m^ = rua 



2^r^ = ZI2Y^ (23) 



p rUa mlLoT 

It is easy to show that for the mean energy Ea = 4.2 keV, axion mass rUa = 17 eV and 
the thickness of the overshoot tachocline zone Lqt = 7.25 • 10^ m the density variations in (23) 
are ~10~^^. It means that the inverse coherent Primakoff effect takes place only when the 
variations of the medium density inside a cylindric volume of the "height" Lqt (see Fig. lib) 
do not exceed the value of ~10~^^. This is a very strong restriction, since it is hard to imagine 
any kind of a physical process in the magnetic flux tube (see Fig. lib) which would "freeze" 
the plasma (low-Z gas) in this magnetic volume so much so that this restriction on the density 
variations is fulfilled. 



^Hereinafter we use rationalized natural units to convert the magnetic field units from Tesla to eV'^, where the 
conversion is 1 T = 195 eV'^ [431. 
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Figure 11: (a) Magnetic loop tubes formation in the tachocline through the shear flows instability development; (b) 
"Capillary" effect in magnetic tubes and the sketch of the axions (red arrovsfs) conversion into 7-quanta inside the 
magnetic flux tubes containing the magnetic steps. Here Lqt is the height of the magnetic shear steps. The tubes' 
rotation is not shown here for the sake of simplicity; (c) Emergence of magnetic flux bundle and coalescence of spots 
to explain the phenomenology of active region emergence (adopted from [59], [69]). 
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In other words, such mechanism that would vahdate the possibihty of such locaUy "frozen" 
plasma existence is not known to us. At the same time, there are some arguments suggesting 
that such limitation is possible. 

First of them is related to the experimentally obtained solar images (Fig. 13, adapted 
from [14]) from the Japanese X-ray telescope Yohkoh (1991-2001) which illustrate the solar 
X-ray activity during the last maximum of the 11-year solar cycle (Fig. 13b). There is currently 
no model alternative to the axion mechanism of sun luminosity which would have described the 
anomalous distribution of the X-ray radiation over the active Sun surface. 

The second argument is related to finding of the axion with mass rua = 17 eV during the 
study of the EBL spectral intensity (see Section 5). On the one hand, this experimentally 
established fact suggests that the solar axions with mass nia = 17 eV are not only a theoretical 
prediction, but they really exist; and on the other hand, it is crucial for substantiation of the 
axion mechanism of Sun luminosity and the solar dynamo ~ geodynamo connection. 

Finally, the third argument is related to the lack of understanding the link between the 
magnetic tubes formation and lifetime in the tachocline, and the length of the solar cycle. It 
means that we do not understand the mechanisms of solar activity as well as the processes of 
generation, accumulation and release of the magnetic energy responsible for the 11-year solar 
cycle as yet. This is applies especially to the current level of understanding of the causes and 
effects of the differential rotation and meridional circulation in the tachocline. Let us remind that 
the tachocline is a thin transitional zone with the width of only O.ObRs, where the latitudinal 
differential rotation of the convective zone turns into almost solid-body rotation of the radiative 
zone (e.g. [101]). It is generally believed that the main process of magnetic field generation — 
solar dynamo - responsible for the 11-year cycle takes place in this zone [102]. However, there 
are no evidence of the 11-year variations in the tachocline so far. Instead, mysterious variations 
of the rotation velocity with period of 1.3 year are observed here [103], and curiously enough, 
they coincide with en estimate of the magnetic flux tube lifetime in the convective zone [68]. 

Although there is no deep and detailed understanding of the magnetic tubes formation in 
the tachocline, one may assume the following picture of this process in the framework of the 
axion mechanism of Sun luminosity. First of all, numerous magnetic tubes appear (Fig. 11a) 
as a consequence of the shear flows instability development in the tachocline. As is shown 
in Appendix B, the pressure inside such tubes is ultralow, which directly leads to formation 
and " floating- up" of the magnetic steps in these tubes (Fig. lib). In other words, a kind of 
"capillary" effect is observed in this case. The whole picture of the magnetic tube spatio- 
temporal evolution in the convective zone depicted at Fig. 11 leaves the question about the 
locally " frozen" plasma existence open, but at the same time it illustrates the process of axions 
conversion into 7-quanta in the overshoot tachocline and thus the mechanism of sun luminosity 
and active solar regions formation in the photosphere (Fig. lie). 

By normalizing the expression (22) the probability of the total conversion of axions into pho- 
tons is assumed to be equal to a unit at given parameters of the magnetic field, i.e. Pa^^y = 1- 
The primary criterion for such choice of the axion-photon coupling strength is the assumption 
about the maximum contribution of the luminosity produced by the axion to 7-quanta con- 
version in the tachocline zone (see Fig. 6c, d) into the total Solar luminosity (Asun) during the 
active phase. 



Aa • ^pr ■ {Ea)pj. + ^Brems ' {Ea) Brems + ^Compt ' {Ea)compt + ^A/1 • {Ea) Ml 

X {AttRIe) X Pa^^ = Asun, (24) 

where^ Aa = 0.90 is a portion of the axion flux that transforms into 7-quanta in the tachocline 
by the inverse Primakoff effect'^; Pa^j = 1; As„„ = 3.84- 10^6 ±1.5% W [86, 106]; Rse = 1.496- 



^In the balance Eq. (24) we considered Ml transition of the ^^Fe nuclei only and ignored the ^^Mn and ^"^Na nuclei, 
because the favorable Boltzman factor of ^^Fe produces the largest cooling rate near Tg ~ 1 [104]. 

''The choice of the Aa value is dictated by the spatial geometry of the solar tachocline zone (see Fig. 6c, d) which 
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10^'^ cm is the distance from the Earth to the Sun; {Ea)pr = 4.2 keV, {Ea)Brems = 1-6 keV , 
{Ea)compt = 5.1 keV, {Ea)Mi = 14.4 keV are the average energies of the solar axions spectra; 
^Pr, ^Brems, ^Compt, ^Mi are the integral solar axions fluxes generated by the Primakoff 
effect, the bremsstrahlung, the Compton process and the Ml transition in the ^^Fe nuclei on 
the Sun respectively, which were obtained using the known differential spectra d^Pr/dE [25], 
d^Brems/dE [6], d^Compt/dE [6], d^Ml/dE [2]: 

$Pr = 3.75 • lO^^c/^^ cm^^g-i^ (25) 

$Brem. = 1.43 • 1035^2^ cm-^ s-\ (26) 

^Corapt = 2.lQ-lQ^^gl, cm-^s-\ (27) 

$Mi = 1.66 • 1023^1 cm-2s-i; (28) 

where ^a^, Qae-, Qan are the axion coupling constants to photons (507)5 electrons [gae) and 
nucleons (ffan) respectively^. 

We obtained only the axion-photon coupling constant ga^y (see (22)) out of the three axion 
coupling constants so far. Next it is not too hard to estimate the axion-nucleon coupling constant 
basing on the hadronic axion models [107, 108]. 

9an = \goP + g?.\, (29) 

where 



rriN 



'' - 2/. 



25 + (3F-Z))l±^^^ 

^ ' l + z + w 

rriN 



[D + F) ^ 



l + z + w 



(30) 



(31) 



Here the value /3 = —1.19 for the Ml transition in the ^^Fe nucleus was calculated in [107, 
108]; spontaneous breaking of the PQ-symmetry (in view of (lO)-(ll)) takes place at the energy 
fa ~ 0.119 . 10^ GeV; mj\[ = 939 MeV is the nucleon mass. The exact values of D and F 
parameters determined from semileptonic hyperon decays are equal to D = 0.808 it 0.006 and 
F = 0.462 lb 0.011 [109]. Parameters z = mu/rrid ~ 0.56 and w = niu/ms ~ 0.029 are quark 
mass ratios [32, 34]. 

The parameter S characterizing the flavor singlet coupling still remains a poorly constrained 
one. Its value varies from S = 0.68 in the naive quark model down to S = —0.09 which is given 
on the basis of the EMC collaboration measurements [110]. The more stringent boundaries 
(0.37 ^ 5 ^ 0.53) and (0.15 ^ 5 ^ 0.5) were found in [111] and [112], accordingly. As a result 
the value of the sum (29) may significantly decrease and, due to negativity of the parameter /3, 
actually vanish. Taking into account that the usually accepted value of u- and d-quark mass 
ratio z = 0.56 can vary in 0.35-^0.6 range [113], the exact interpretation of experimental results 
is significantly restricted. 

Calculations performed using the expressions (29)-(31) and (13) let one derive the value of 
gan = 2.46 . 10~^ for the axion-nucleon coupling constant at 5 = 0.55. However, it should be 
noted that such a model-dependent value does not suit our needs, and here is why. 



is described in more detail below (see Fig. 9). Let us note that the portion of the "invisible" axions equal to (1 — Aq) 
satisfies the neutrino limit on axions, i.e. the Gondolo-RafFelt criterion [82, 105] at the same time. 

^Although Eqs. (25)-(28) were obtained in the framework of the standard model of the Sun, they may be applied 
to the axion model as well, since the basic thermodynamic parameters (temperature, pressure, plasma density etc.) of 
the solar core are roughly the same in both models, and almost all the axions are produced in the solar core, regardless 
of the exact way of their birth. 
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It is easy to show using (28) that the luminosity Aj\,f i of the axions produced by Ml transition 
of ^^Fe nuclei depends on the solar photon luminosity A^^^ in the following way^: 

Ami = 2.8 • lO^^LA^.n, (32) 

Substituting the value of gan = 2.46 • 10~^ into (32) we derive that the relative axion 
luminosity Ami is about ~ 2%. This value is inadmissible for the axion mechanism of Sun 
luminosity, because otherwise the resulting solar photon spectrum (Fig. 7) would contain a 
rather high peak near Ea = 14.4 keV. As a matter of fact, if it is there, then it is very weak 
considering the spectrum uncertainty in this band (see Fig. 9 in [115]). 

In this connection let us from now on assume that the relative axion luminosity Ajv/i is 

Aa/i/A5„„ = 0.003. (33) 

The choice of such relation was made so that the relative axion luminosity Ami allowed 
the 14.4 keV peak existence within the experimental Sun photon spectrum uncertainty and 
conformed to the axion-mucleon coupling constant from the theoretically allowed limitations 
known as the SN1987A limit (3 • lO"'^ ^ gan ^ lO'^ [116, 117]). 

Substituting (33) into (32) we obtain a consistent value of the axion-nucleon coupling con- 
stant i° 

gan = 3.2 • 10-^ (34) 

Now let us turn again to Eq. (22) which describes the axion mechanism of Sun luminosity hy- 
pothesis. Apparently, the found axion coupling constants to photons {ga-y = 7.07-10^^^ GeV^^) 
and nucleons {gan = 3.2 • 10~^) let one calculate the axion-electron coupling constant (gae) by 
means of Eqs. (22)-(28). It turned out to be 

gae = 5.28 -IQ-^K (35) 

It is interesting to note that within such approach, bremsstrahlung (~ 2.59 • 10^^ W) has the 
major part in the Solar luminosity, which is ~ 67.48%, while the Compton process (~ 1.24- 10^^ 
W), the Primakoff effect (~ 3.19 • 10^3 W) and the Ml transition of "Fe nuclei (~ 9.93 • 10^2 
W) make 32.41%, 0.08% and 0.03% respectively. Differential solar axions spectra as observed 
at the Earth originating from bremsstrahlung, the Compton process, the Primakoff effect and 
Ml ground-state nuclear transition in solar ^ Fe are shown on Fig. 2b. 

Thus, it is obvious that the mechanism of luminosity and the X-ray spectrum shape for the 
active and quiet Sun (Fig. 7) can be easily explained by our model, in which axions are converted 
into 7-quanta in the solar tachocline under certain conditions. First of all, the 7-quanta energy 
spectrum generated by axions in the solar tachocline zone due to the channeling effect practically 
does not change to the boundary of the photosphere of the Sun, where it transforms because 
of the Compton scattering, as is shown in Figs. 6c, d and Fig. 7. And, finally, it is obvious that 
the integral of this spectrum (see Fig. 7) by virtue of the equality (24) coincides by the order 
of magnitude with the estimation of the Sun luminosity. 

2.5 Invisible axions and Solar Equator effect 

As it was already stated before, the solar magnetic field variations (Fig. 3) must drive (with 
the help of the solar axions and the inverse Primakoff effect) the total solar irradiance (TSI) 
and 14.4 keV axions variations at the same time (see the inset in Fig. 12). It is important to 



^Here we follow the calculations by Derbin (e.g. [3]), that is why the expression (32) differs slightly from the 
analogous expression (2.13) in [114]. The difference mainly comes from the different values of the Doppler spectrum 
broadening used during the calculation of the monoenergetic solar axions flux produced by the ^"^Fe nuclear de- 
excitations (e.g. [3, 114]). 

^"^Hereinafter by gan and gae we always mean \gan\ and \gae\ respectively. 
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Figure 12: Time evolution of Pa^j and TSI during 1975-2010. Inset: time evolution of (a) the variations of the 
magnetic flux in the tachocline zone of the Sun [12]), (b) TSI annual variations [106]. Curves are smoothed by the 
sliding intervals in 5 and 11 years. 

note here that, however strange it may seem, TSI variations are not the modulator of the Earth 
climatic system (ECS) global temperature, because the strong inverse^^ correlation with the 
22-year lag is observed between them [9]. And vice versa, the variations of the solar axion flux 
manifest the strong positive correlation with the global temperature variations with the same 
time lag. This fact plays a key role in the new global climate theory [118-120], which considers 
the variations of the 14.4 keV solar axions (which are resonantly absorbed in the Earth core) a 
trigger-like modulator of all the thermal processes in ECS and, particularly, in the atmosphere. 
However, this problem requires a special discussion and therefore will be examined in a separate 
paper. 

On the other hand, as follows from Fig. 12, the TSI variations during the active phase of 
the Sun are so small (~ 1 W/m? [106]), that the relative portion of 14.4 keV axions must also 
be small at the Earth. 



PTSI 



ARa _ (TSI) ■ 47rr|^ 



Ra ^Sun 

Therefore their heat power in the Earth core 



10' 



(36) 



T^ = PTSI ■ Ra 



A^Fe • Ea 



20 VF, 



(37) 



is not enough not only for the geomagnetic field generation (which requires at least ^0.1 TW [18]), 
but for the geomagnetic field variations also (~0.01 TW). Here Ra = 5.16-10~^5^„, NpJ^ ~ 3-10 



l47 



is the number of ^^Fe nuclei in the Earth core, Ea = 14.4 keV is the ^^Fe solar axions energy. 

^^A lot of climatologists still believe that it is the TSI variations that are responsible for the Earth global temperature 
variations obstinately disregarding the facts of a very small contribution of TSI variations into the energy balance in 
the Earth atmosphere and the inverse correlation between TSI and global temperature variations (see Fig.l in [9], 
where the ocean level is a proxy for the global temperature). 
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Figure 13: Top: Solar images at photon energies from 250 eV up to a few keV from the Japanese X-ray telescope 
Yohkoh (1991-2001) (adapted from [14]). The following is shown: on the left, a composite of 49 of the quietest solar 
periods during the solar minimum in 1996. On the right, the solar X-ray activity during the last maximum of the 
11-year solar cycle. According to Fig. 6 (and Fig. 13, bottom) most of the X-ray solar activity (right) occurs at a 
wide bandwidth of ±45° in latitude, being homogeneous in longitude. Note that ^ 95% of the solar magnetic activity 
covers this bandwidth [121] (see also a similar topology for microflares measured with RHESSI [122]). 
Bottom: Schematic picture of the solar tachocline zone, the Earth's liquid outer (the red region) and inner (the brown 
region) core. Blue lines on the Sun designate the magnetic field. In the tachocline axions are converted into 7-quanta, 
which form the experimentally observed solar photon spectrum (Fig. 7) after passing through the photosphere. Solar 
axions moving towards the poles (blue cones) and in the equatorial plane (blue bandwidth) are not transformed by the 
Primakoff effect, since the magnetic field vector is almost coUinear to their momentum vector in these regions. Solar 
axions are then resonantly absorbed by iron in the Earth core transforming into 7-quanta, which are the supplementary 
energy source in the Earth core (see the text). 

At the same time it is not difficult to see that the relative part of the axions A^ that are 
almost not affected by the Primakoff effect in the polar (Apoi) and equatorial (Aggu) sectors of 
the tachocline zone (Fig. 13) is a considerable quantity: 

A„ = A°^„ + A^,„ (38) 

If we assume that the equatorial surface (Sequ) formed by two cones going from the center 
of the Sun, and a part of the sphere with the radius of the Sun (RSun) has the dihedral angle 
of ~ 5°, it is easy to show that^^ 



^This estimate was made on the basis of the numerous computational experiments on magnetic field evolution in 
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K,u = 1^ ~ 0.05. (39) 

The expression (39) means that because of the quasi-colhnearity of the Earth and Sun 
rotation axes, the main axion flux directed towards the Earth originates from the equatorial 
sector of the Sun (Fig. 13). In this connection a short remark should be made regarding 
the anti-correlation between the solar magnetic field and geomagnetic field variations. Weak 
variations of TSI are obviously produced by the solar magnetic field variations. The same 
solar magnetic field variations are the cause of the "equatorial" axion flux variations. In other 
words, the "equatorial" effect not only generates the "invisible" axions, but also modulates 
their intensity inversely proportional to the solar magnetic field changes, producing the observed 
inverse correlation between them (Fig. 11). The latter is supposed to be the main cause of the 
anticorrelation between the solar magnetic field variations and the geomagnetic field variations. 

The assumptions used to derive the expressions (16) and (39) are called forth by the necessity 
of justification both the axion mechanism of Sun luminosity and the axion mechanism of solar 
dynamo - geodynamo connection, and will be additionally substantiated below. 

2.6 Power required to maintain the Earth magnetic field and 
nuclear georeactor 

It is not hard to show that the resonant absorption rate of 14.4 keV solar axions in the Earth 
core, which contains the N^!^ nuclei of ^^Fe isotope, is about [3] 

Ra = h.lQ-l^-^gl^-N%-A%^, (40) 

where Ae^^ is the portion of axions reaching the Earth via the solar equator effect (see (39)). 

It is known, that the number of ^"^Fe nuclei in the Earth core is Np^ ~ 3 • 10^^ and the 
average energy of ^^Fe solar axions is Ea = 14.4 keV. Then with an allowance for Eq. (40) and 
the value of the axion-nucleon coupling constant (35) the maximum energy release rate AD"^^^^ 
in the Earth core is equal to 

AD„\„,, = i?,.K^1.2 kW. (41) 

This estimate of the heat power supplied to the Earth core by the absorbed axions (41) 
apparently is much less then the value necessary to generate the magnetic field of the Earth 
(^ 0.1 TW [18]). Moreover, it is small even in comparison with the heat power fluctuations 
(~0.01 TW) responsible for the geomagnetic field variations in the Earth core (see the inset in 
Fig. 12). 

It is easy to illustrate this using the known dependence of the core magnetic field Bin on 
ohmic dissipation Dohmic in the Earth core [18, 84, 123] in the form: 



Dohmic = ^^Bt^Q.l TW, (42) 

where t/ ~ 2 m? / s is magnetic diffusivity [17], V = (4/3)7r-R^ is the volume of the Earth core, 
Re = 3480 km is the radius of the Earth core, /u ~ 1 is permeability, Ib ~ 0.8 • 10^ m is 
the characteristic length scale on which the field vector changes [17], Bin ~ 4 mT is the core 
magnetic field [124]. 

In order to estimate the heat power fiuctuations ADohmic responsible for the magnetic field 
variations in the Earth core let us represent (42) in differential form: 

2ry • V 

U-J-J ohmic ^ P}~ Oindtiin, V^'^) 



the convective zone of the Sun (e.g. [12] and Refs. therein). 
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The right-hand side of (43) contains the known estimates except for the magnetic field 
fluctuations Ai3j„, ~ dBin- On the other hand, there are well known long-term records of 
magnetic field variations measured on the Earth surface (e.g. [13]). This lets us estimate dDo^^ic 
by writing down the expression (40) in the following form: 

<;Z>„.„„ = ^B.„=§i, ,44) 

jjL at 

where the magnetic field variations (dBin/dt) in the Earth core appear after taking into account 
the magnetic diffusion 

dB,r,^^5t = ^Jl. (45) 

dt dt 1] ^ ^ 

If one also takes into account the known relation between the inner (dBin/dt) and outer 
[dBout/dt) geomagnetic field variations [125], 

then it is possible to make an estimate of the ohmic dissipation fluctuations {IS.Dohrnic ~ 
dDohmic) basing on (43)- (45) necessary for inducing a certain number N of Taylor cells in 
the core [125]. 



ADoWc = —B^n ( ^ ) '^^ ^ 0.02 TW, (47) 



2y„ / Re Y dBout 
H \jRcJ dt 



Here Re = 6357 km is the radius of the Earth, 7 = 0.1 [125], {dBout/dt) ~ 20 nT/yr is the 
annual variation of the external magnetic field of the Earth core [20, 126]. 

A natural question arises from the stated above about the way that 14.4 keV solar axions 
may provide an effective mechanism of solar dynamo - geodynamo connection while supplying 
a rather low heat power. In other words, how does this problem reduce to the mechanism 
of small heat perturbations critical influence on the convective process in the Earth liquid 
core. The problem is stated this way because if there is an effective mechanism of convective 
instabilities generation by weak heat perturbations in the Earth liquid core, then this effect 
may simultaneously cause substantial weakening of the convective heat removal from the Earth 
solid core. The intense weakening of the heat removal from the Earth solid core surface, in its 
turn, causes the corresponding temperature increase in the solid core boundary layer. This is 
very important because it promotes the subsequent effective convection stability recovery in the 
liquid core. 

There are strong grounds to believe that there is a natural nuclear georeactor operating at 
the boundary between the Earth solid core and liquid core. The analysis of the KamLAND 
experiment neutrino spectra for 2002-2008 shows that the heat power of this non-stationary 
traveling wave reactor (TWR) is about 30 TW [127-129] . The heat power of such TWR depends 
on the nuclear fuel composition and the medium temperature. The latter is because of the fact 
that according to [128, 129], 238-jj g^j^^ 239p^ capture and fission cross-sections depend quasi- 
linearly on the temperature of the neutron-multiplicating medium in the 3000-5500 K range 
(Fig. 14). 

These peculiarities of TWR are responsible for the positive feedback that leads to the re- 
actor heat power increase after the corresponding boundary layer temperature increase (see 
Fig. 14). The georeactor heat power growth lasts until the steady heat removal from the TWR 
is reestablished, which implies restoration and stabilization of the convection in the Earth liquid 
core. 

Now let us turn back to the physical essence of the axion mechanism of the weak convective 
instability thermal perturbations in the liquid core. In this connection it should be noted that 
the convection in the Earth core is compositional. It means that there are some light elements 
originating, in particular, from the ^^^Pu nuclei fission which take part in the convection along 
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Figure 14: Dependence of (a) capture cross-sections and (b) fission cross-sections for ^35^ (bi^j^e)^ 238^ (^gY&ea), and 
239p^ (red) averaged over the neutron spectrum on the fuel medium temperature for the limiting energy (3kT) of the 
Fermi and Maxwell spectra [129]. The neutron spectra averaging procedure was applied for the concentrational fuel 
composition of the nuclear georeactor discussed in [127]. 

with the "iron" component [127-129]. It turns out that the convective instability may appear 
in such media even under hydrostatically stable density stratification, i.e. when the density 
decreases with height [130-132]. 

It is known that the phenomenon of the convective instabiHty caused by the double (differ- 
ential) diffusion was discovered rather long ago and has been described in numerous overviews 
and monographs in detail (e.g. [130-132]). The principal role in this case usually belongs to 
the difference between the two hydrodynamic components of heat and admixture [130]. The 
convection caused by double diffusion is generally believed to appear when the thermal medium 
stratification is stable, while the weakly diffusing admixture (e.g. light elements) introduces 
a destabilizing contribution into the density stratification. Although this contribution may be 
relatively small, it may be enough for destabilization of a stably stratified (in terms of density) 
system owing to the mentioned effects [132]. 

However, we are interested in the conditions of the convective instability formation in a qual- 
itatively different situation, particularly, when a weakly diffusing admixture, on the contrary, 
introduces a stabilizing contribution into the density stratification. This contribution may even 
exceed the thermal instability in absolute value. Such possibility may seem paradoxical at first 
glance since due to double diffusion effects the slowly diffusing admixture usually has the much 
greater impact on the convective instability than the quickly transportable heat, all other fac- 
tors being equal. Let us show that this is not always the case by analyzing the situation when 
there is a slow background motion along the gravity force described in [132] in detail. 

A problem on convection on the background of the slow (relative to the characteristic speed 
of the studied convective motions) vertical motion was first considered in [133]. It is of a con- 
siderable interest for us, since the resonant absorption of 14.4 keV solar axions in the iron nuclei 
may be considered as a process that induces a slow descending background motion (Fig. 15a) 
in the convective medium of the liquid core. 

Following [132], let us consider a single-component medium with its density depending on 
the temperature T only (neglecting the admixture stratification effects). In other words, we are 
considering a modification of the classical Rayleigh-Benard problem on the convective stability 
of a liquid between two horizontal plates [130]. A slow vertical motion along the gravity force 
is assumed to be present in the background state. For the sake of simplicity let us consider 
a motion with the velocity VF < independent of the vertical coordinate z counted from the 
top boundary. Let us also consider the bottom and top boundaries temperatures T},ot and Ttop 
fixed and denote the difference between them by AT. Heat transfer in the background flow is 
described by the equation 
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Figure 15: (a) A sketch of the thermal convection with the resonant absortion of 14.4 keV solar axions by iron nuclei 
producing the 14.4 keV 7-quanta flux (red arrows). Here 7-quanta flux emulates a slow background descending motion 
in the convecting medium. The green arrow denotes the convection direction, (b) Distortion of the vertical background 
temperature distribution by the downward motion: (1) the zero background vertical velocity, w = H/h = 0; (2) w = 10 
(adapted from [132]). 
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d?T 



where k is the thermal diffusivity. A solution for the two boundary conditions mentioned above 
may be written down in the form 



G(z) 



exp(— ^) — exp(— w) 



1 — exp(- 



-w] 



(49) 



Here 0(z) = \T{z) — Ttop]/ ^T is the dimensionless temperature deviation, ^ = h/ k is the 
dimensionless vertical coordinate, and h = k/W is the reference height associated with the 
vertical motion (infinity in the quiescent fluid). The key dimensional parameter is 



w 



H 



W{H/k) 



(50) 



where H is the fluid layer thickness. In the absence of the background vertical motion (in the 
limit W ^ 0, h -^ 00, and w; — >■ 0), the result is the expected linear proflle 



9 = 1- z/H, 



(51) 



i.e. the solution whose stability is analyzed in the classical Rayleigh-Benard problem. Fig. 15b 
shows the vertical proflles of (51) for w = and w = 10 [132, 133]. Apparently, the background 
medium sinking "pushes" almost all temperature difference AT to the bottom boundary where 
it concentrates within a layer with the thickness h = kJW . 

Generally speaking, a rigorous stability study of the stationary state with a curvilinear 
temperature proflle and the background sinking is a rather complex problem. An estimate 
of the effective Rayleigh number for the bottom sublayer, which incorporates virtually all the 
vertical temperature difference AT (Fig. 15b), was made on the basis of the simple physical 
reasoning in the paper [133]. 



Ra 



a 



g- AT -h^ a-g- AT 



Here a is the thermal expansion coefflcient of the fluid, i^ is the kinematic viscosity, and g 
is the gravitational acceleration. 
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Let us denote the value of the effective Rayleigh number, which corresponds to the loss of 
stability, by Racr- The sinking rate sufficient for stability loss prevention in this case is expressed 
by the equation 

W^r-( '"''f"'' X^'-W-' m.s-\ (53) 



\ U ■ Racr / 

where Rucr is a complex function of e [134] 

Racr ~ E^-^^ 0.21e~2 + 22.4(1 - e^)^/^ ~ 10"^ e = Rir,/Rout, (54) 

For example, setting k ~ 0.1 m? • s, z/ ~ 10'^ m? ■ s (effective turbulent transport coefficients 
characteristic for the liquid core [135]), a ~ 10"^ K'^ [136], AT ~ 1000 K [136], 5 = 11 m • 
s~^ [136], E ~ 10^ is the Ekman number for a given k [135], one obtains Wcr ^ 10~^ m ■ s~^. 

It is interesting to note that the background sinking leads to an effective Rayleigh number 
decrease, i.e. Racr < Rci, where Ra > 10^ is Rayleigh number in Earth core [136], and conse- 
quently, to a decrease in convective instability formation probability. This value of downward 
velocity is also in good agreement with the value of the characteristic velocity of the Earth 
core convection (~ 4 • 10~^ m • s~^ [135]), while the results obtained in both field experiments 
and numerical simulations (see Fig. 15b) demonstrate that downward flow with this velocity 
suppresses convection [133]. 

Now let us pass on to a two-component medium which has its unstable thermal stratification 
under the absence of the vertical motion overcompensated by a stable admixture stratification. 
As is shown in [133], if there is a slowly diffusing admixture stratification along with the tem- 
perature stratification (light elements in our case), admixture is suppressed by the background 
movement more effectively than the heat in such two-component medium. In other words, even 
the presence of a very slow downward movement may pull the admixture down as opposed to 
the heat and thus cancel its stabilizing effect, and the system becomes unstable. It is also known 
that the temperature profile may be deformed as well under more intense background motions. 
It is important to keep in mind that the mentioned effects are possible under very small vertical 
velocities of the background motions. As the authors of [133] point out, we are dealing with a 
new kind of instability. Note, however, that the times required for a system to evolve into the 
unstable steady states considered above may be very large [133]. 

It may be concluded that the vertical background motions may prevent the convective in- 
stability formation in the single-component media while leading to a destabilization of a two- 
component medium layer. It happens because the background motions are an immediate cause 
of the effective vertical drift of a slowly diffusing admixture. An important fact to remember is 
that the above-mentioned effects are possible even under very small vertical background motion 
velocities''^. 

The problem on convective instability caused by the background motion is examined in 
its simplest form so far. The rotation effects, magnetic field influence, nonlinear background 
motion velocity etc. were not taken into account here, but all of them are actually present in a 
traditional composite media magnetohydrodynamics in the Earth core. A detailed consideration 
of these effects is beyond the scope of the present paper. The purpose of the current section 
is to demonstrate a possibility of a nontrivial impact of background motions, which may be 
produced by the resonant absorption of the "iron" solar axions along with the convection in the 
Earth liquid core, within a simple model. 



^^It is interesting to note here that a problem on convection in presence of the slow background motions is extremely 
urgent for the known geophysical applications related to e.g. cloud patterns and atmospheric circulation [133, 137- 
139]. For example, the atmospheric and oceanic convection often takes place on the background of the processes 
with much larger horizontal scales (cyclones and anticyclones) which are characterized by the average vertical motions 
several orders of magnitude slower than those that appear during a convective instability formation. According to the 
natural experiments [139], even a slow background sinking of the medium can effectively suppress the convection in 
the atmosphere. 
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Thus, the essence of the axion mechanism of solar dynamo - geodynamo connection hes in 
the fohowing. The resonant absorption of 14.4 keV solar axions by the iron of the Earth core 
induces a vertical background motion along the gravity force (Fig. 15a), which in its turn "pulls" 
almost all the temperature difference AT down to the bottom of the liquid core (Fig. 15b) and 
concentrates it within a layer of the thickness h = k/W. As it was noted above, this effect takes 
place both in single-component and in two-component media. 

An important result of these processes is a substantial attenuation of a heat removal from the 
Earth solid core surface which leads to a corresponding temperature increase in the boundary 
layer between the liquid core and the solid core where the nuclear georeactor (TWR) resides. 
As it was shown in [129], one of the peculiarities of such TWR is that its heat power output 
depends both on the fuel composition and the medium temperature. It means that increase 
of the temperature in the boundary layer at the solid core and liquid core border leads to 
a corresponding increase of the nuclear georeactor power output (see Fig. 14). As a result, 
the georeactor heat power output grows until a steady heat removal is re-established, i.e. the 
convection is re-established and stabilized in the liquid core (up to the "next" variation of the 
thermal perturbations by axions!). 

Therefore if such axion mechanism of solar dynamo - geodynamo connection exists, then 
the ohmic dissipation caused by a resonant 14.4 keV solar axions absorption in the Earth core 
should be connected with the heat power perturbations ADohmic, responsible for the magnetic 
field variations in the earth core, by the following relations: 

Bin = ^- Bf^, (55) 

dUohmic = ? "-L^o/imio (^o) 

where the trigger gain ^ in our case (see (38) and (44)) is equal to 

^~ 3.4 -10^. (57) 

Here Bin and Sf^ are the magnetic fields in the Earth liquid core produced by the nuclear 
georeactor and the solar axions respectively. The physical sense of the expressions (55)- (56) 
reveals the reason why all of the known candidates for an energy source of the Earth magnetic 
field [15] cannot in principle explain one of the most remarkable phenomena in solar-terrestrial 
physics - a strong (inverse) correlation between the temporal variations of magnetic flux in the 
overshoot tachocline zone [12] and the Earth magnetic field (Y-component) [13] (Fig. 3). 

3 Axion mechanism of Sun luminosity and CUORE 
experiment 

Recently a CUORE-experimental search was performed for axions from the solar core from 
14.4 keV Ml ground-state nuclear transition in ^^Fe [1]. The detection technique employed a 
search for a peak in the energy spectrum at 14.4 keV when the axion is absorbed by an electron 
via the axio-electric effect. The cross section for this process is proportional to the photo-electric 
absorption cross section for photons [140]: 

87TaEM \ fa J yrUeC^ J 

where Xg Ri 1 , mgC^ is the electron mass in GeV, uem = 1/137, and Upe = 55.336 cm? ■ gm~^ 
is the photoelectric cross sections for Te02 (taken from [141]). 

Substituting the values of these constants one may rewrite Eq. (58) in the following form: 

2.18-10-11 GeV^ f BaV ,,„, 

aae = J, ^^J ape, (59) 
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where Ea is in keV and fa is the Peccei-Quinn scale in GeV. 

Hence the estimated absorption rate of 14.4 keV solar axions N quo RE detected by the 
axio-electric effect (58) in Te02-detector [m = 3 kg) of the CUORE experiment is 

^a-T^CTae-eo ^ NcuoRE = ^-^^ count ■ kg^^d^^ , (60) 

where ed is the detection efficiencies; the flux ^a at the Earth (28) is 

$a = 1.66 • 1023(5,„)2 cm-h-\ (61) 

It is the Eq. (60) that made it possible for CUORE collaboration to place a bound (at 
S = 0.55) on the axion coupling constant of /« ^ 0.76 • 10^ GeV at 95% C.L. (Fig. 16b). 
According to Eq. (10), the limit on fa translates into a mass limit nia < 8 eV. 

It is necessary to note that if one takes into account the axion mechanism of Sun luminosity 
and solar dynamo-geodynamo connection, expression (60) with respect to solar equator (A" r^ 
0.05), becomes 

^ . NcuORE 0.63 -1^-1 /co^ 

$a • vTCTae • Ed ^ —T^ = "T^— counts ■ kg d , (62) 

^equ ^equ 

In other words, within the framework of such mechanism, it is necessary to keep in mind 
that because of the solar equator effect, only a part of the total axion flux $a (61) equal to 
^equ^a arrives to the Earth. The solution (62) in this case lets us use the expressions^^ [1] 

o /6.2-106 Gey\ A^ ^A-z\ 

93 = -7.8. 10- ( — ^^ — j((^+^)rT^)- ("> 

to determine (see Fig. 16a) the /* for a value of S" = 0.55^^ 

/* ^ 0.353 • 10^ GeV, (65) 

calculate (see Fig. 16a) the value of axion-nucleon coupling constant (at S = 0.55) 

51 ^ 4.8 • 10-^ (66) 

and estimate the axion mass by means of (10) and (65): 

ml ^ 17 eV. (67) 

The estimates (66) and (67) demonstrate and excellent agreement with the axion-nucleon 
coupling constant (34) and axion mass (11) obtained in the framework of axion mechanism of 
Sun luminosity and solar dynamo - geodynamo connection. 

These results make it possible to estimate the limit on the axion-electron coupling constant 
gae as well. Derevianko et al. [143] and Derbin et al. [7] showed that the cross section aae{Ea) 
for the axio-electric effect is proportional to the photoelectric cross section Upe {E) , and is given 
by the formula [7]: 

.„(E)=.,(E)f3A_(j?^y(i_0, (08) 

^^Let us point out that that the expressions (63) and (64) used in CUORE experiment data processing are slightly 
different from the corresponding expressions (30)-(31) derived in [107, 108]. 

^^It is appropriate to mention here that E143-experinient provided the value for the parameter S equal to 0.30 ± 
0.06 [142]. 
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Figure 16: Expected rate in the axion region as a function of the fa axion constant for different values of the nuclear 
S parameter. The horizontal line indicates the upper limit obtained in the present paper (fa ^ 0.353 • 10^ GeV for 
S — 0.55) and in CUORE-experiment {fa ^ 0.76 • 10^ GeV for S = 0.55) [IJ.The insets show the correspondence 
between the axion-nucleon coupling constant and Peccei-Quinn energy scale (a) in our model and (b) in CUORE- 
experiment. Bounds are given for the interval 0.15 ^ S* ^ 0.55. 

where Ea is the axion total energy, f3 is the axion velocity divided by the velocity of light, and 
Qae is the dimensionless axion-electron coupling constant. At /? — )• 1 and /? — )• 0, this formula 
coincides with the cross sections for relativistic and nonrelativistic axions obtained in [144]. 

Since the axion mass is small (see (67)), let us consider a relativistic form of Eq. (68) (i.e. the 
case of /? — 7- 1) and compare it to the analogous expression (58). As a result of such comparison 
we derive the upper limit on axion-electron coupling constant with respect to (65): 



2x'„mp 

< 2.89 



f* 

J a 



10" 



(69) 



The limit (69) obviously does not contradict the value of axion-electron coupling con- 
stant (35) obtained in the framework of axion mechanism of Sun luminosity. 

So it may be concluded that the new estimations for the strength of the axion-photon 
coupling {ga-y ~ 7.07 • 10~^^ GeV~^), the axion-nucleon coupling {gan ~ 3.2 • 10""^), the axion- 
electron coupling {gae ~ 5.28-10~^^) and the axion mass {ma ~ 17 eV) obtained in the framework 
of axion mechanism of Sun luminosity and solar dynamo - geodynamo connection are in good 
agreement with the CUORE experiment data. It brings hope that this hypothesis will be 
justified by the future CUORE experiment with the expected exposure of 1.4 • 10^ kg ■ day which 
is significantly larger than the current one (43.65 kg ■ day) [1]. 

4 Axion mechanism of Sun luminosity and other 
important experiments 

In view of the axion mechanism of Sun luminosity, let us analyze the data of known experiments 
on measuring the axion coupling to photon, nucleon and electron for different axion mass ranges. 

4.1 Axion coupling to a photon 

Fig. 17a shows virtually all the major experiments on estimating the limits of the axion-photon 
coupling constant. The statement of the problem in all presented experiments included the 
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Figure 17: (a) Exclusion regions in the "m^ - 5a^"-plane achieved by CAST in the vacuum [26, 145], ^He 
and ^He phase [146]. We also show constraints from the Tokyo hehoscope [147-149], BNL telescope [150], SO- 
LAX [151], COSME [152], CDMS [153], DAMA [154] and the hot dark matter limit (HDM) for hadronic axions 
TOa < 1-05 eV jc? [155] inferred from WMAP observations of the cosmological large-scale structure. The yellow band 
represents typical theoretical models with \ElN — 1.95] = 0.07 -^ 7. The red solid line corresponds to EjN = (the 
KSVZ model). The field theoretic expectations are shown together with the string theory Z12-1 model of Choi, Kim, 
and Kim (the green line) [156, 157]. (b) Same as (a), but all the experimental data are corrected with account for the 
solar equator effect and a contribution of the Primakoff effect into the total Sun luminosity. A red star denotes the 
result obtained in the present paper. 

measurement of the axion flux that could be produced in the Sun by the Primakoff conversion 
of the thermal photons in the electric and magnetic fields of the solar plasma. The difference 
between these experiments consisted in the axion flux detection technique only: by the axion- 
to-photon reconversion (the inverse Primakoff effect) in laboratory transverse magnetic [26, 114, 
145-150] and electric (the intense Coulomb field of nuclei in a crystal lattice of the detector plus 
the Bragg scattering technique [98-100, 104, 151, 152]) fields. 

The new constants for the axion mechanism of Sun luminosity, obviously, should be calcu- 
lated with due regard for the solar equator effect and Primakoff effect contribution into the total 
solar luminosity. It is not hard to show that in this case they should be 



(9. 



a-yJ 



^Brems + ^Compt + ^Pr + ^AIl 



$ 



Pr 



equ 



9a7 



(70) 



Hence, the values of the partial components of Sun luminosity {^Brems ~ 67.48%, ^compt ~ 
32.41%, ^Pr ~ 0.08%, and $mi ~ 0.03%) and the solar equator effect probability {Aggu ~ 0.05) 
lead to a new value of g* constant: 



g* = 1.26- ga-y. 



(71) 



Fig. 17b shows the new limits obtained for the axion mechanism of Sun luminosity taking 
into account (71) for the corresponding experiments. The change relative to the initial picture 
(Fig. 17a) is rather small, since the small contribution of the Primakoff effect into total Sun 
luminosity "effectively" compensates the influence of the solar equator effect. At the same 
time it is interesting to note that our values for the axion-photon coupling constant (5*^ ~ 
7.07 • lO"^"^ GeV~^) and the axion mass (~17 eV) fit well into the existing limits (a red star in 
Fig.l7). 
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Figure 18: (a) Summary of limits on axion-electron coupling. The limits shown include the astrophysical bound from 
the solar neutrino flux [82], dedicated axion experiments by Derbin using '^^^Tm [6] and Si(Li) [7]; CDMS, CoGeNT, 
DAMA and XMASS data obtained from [153, 158-160]; (b) the data from (a) with a due correction associated with 
the solar equator effect and a contribution of bremsstrahlung and the Compton process into the total Sun luminosity, 
except for the astrophysical bound from the solar neutrino flux which is corrected for the contribution of the Compton 
process into total luminosity only. 



4.2 Axion coupling to an electron 

Let us now consider the paper by Derbin et al. [7]. In this paper the axio-electric effect in 
siUcon atoms is sought for solar axions appearing owing to bremsstrahlung and the Compton 
process. Axions are detected using a Si(Li) detector placed in a low-background setup. As a 
result, new model-independent constraints have been obtained for the axion-electron coupling 
constant and the mass of the axion. For axions with a mass smaller than 1 keV, the resulting 
bound is Qae ^ 2.2 • 10^^° (at 90% C.L.). 

Obviously, the new constants for the axion mechanism of Sun luminosity should be derived 
with due account taken of the solar equator effect and a contribution of bremsstrahlung and the 
Compton process into the total Sun luminosity. It is not hard to show that in this case they 
have the following form: 
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(72) 



irems ~r ^Compt 

This time, the values of the partial components of Sun luminosity {^Brems ~ 67.48%, 
^Compt ~ 32.41%, ^Pr ~ 0.08%, and <I>mi ~ 0.03%) and the solar equator effect probabil- 
ity {Aequ ~ 0.05) lead to a new estimate of the g*f, constant limit for the axions with a mass 
smaller than 1 keV: 



9*ae ^ 4.47 • ga 



9.8-10 



"10 



(73) 



where gae ^ 2.2 • lO^^^ [7]. 

Fig. 18b shows the new limitations obtained in the context of the axion mechanism of 
Sun luminosity for the corresponding experiments, according to (72). There is a substantial 
change in Fig. 18b relative to the initial Fig. 18a, because the relatively high contribution of 
bremsstrahlung and the Compton process into the total Sun luminosity cannot compensate the 
solar equator effect probability largely. At the same time, the determined values of the axion- 
electron coupling constant {gae ~ 5.28 • 10~^^) and the axion mass (~ 17 eV) fit rather well into 
the known limits (a red star in Fig. 18) 
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5 Axion dark matter and extragalactic background 
light 

Naturally, having the complete axion "portrait", we are entitled to ask a question about the 
conditions of its detectability. In other words, provided that axions do exist, once they are 
produced, where are they likely to be found? Let us turn to the astrophysical observations. 
Since the current temperature of the cosmic microwave background radiation is T = 2.35 • 
10^'* eV [161], axions are nonrelativistic and have been since before decoupling. Therefore axions 
should, in accord with the equivalence principle, fall with baryons and any other particles into 
the various potential wells which develop in the Universe [162] . The most likely place to find light 
relic axions is in clusters of galaxies and the halos of galaxies. It is always possible, however, 
that the lines of sight to localized regions are partially obscured by previously unrecognized 
amount of absorbing material, so that the decay photon flux is underestimated [163]. To get 
around this problem, we shall consider the diffuse extragalactic background (EBL) rather than 
the flux from any particular region of the sky. 

Before we proceed to the discussion of the axions emission and detectability questions, let 
us consider some important theoretical limitations on the axion parameters. So, within the 
framework of our mechanism, the new estimations of the strength of the axion coupling to 
a photon (^^^ ~ 7.07 • 10~^^ GeV~^), the axion-nucleon coupling {gan ~ 3.2 • 10~^) and the 
axion-electron coupling {gae ~ 5.28 • 10^^^) have been obtained. It is necessary to note that 
obtained estimations cannot be excluded by the existing experimental data (see Fig. 17 and 
Fig. 18), because the discussed above effect of solar axion intensity modulation by magnetic 
field variations in the solar tachocline zone was not taken into account in these observations. 
The obtained estimates for the strength of the axion-photon coupling and the axion-nucleon 
coupling also cannot be ruled out by the existing theoretical limitations known as the globular 
cluster star limit {gaj < 6 • 10~^^ GeV~^) and the red giant star limit {gae < 3 • 10~^^) [105], 
since these values are highly model-dependent^^. That is to say these limitations have a high 
level of uncertainty because of the absence of the standard theoretical model for globular cluster 
stars and red giant stars (see, for example the analysis of theoretical models of the hot core 
in [105] the cross section for axion absorption in [117], and a note for Fig. 3 in [105]). 

There is one more important fact related to the so-called "axion trapping" effect which was 
not taken into account quantitatively in the paper by Raffelt [105]. It is known [164-169], that 
the axion flux from the supernova can be suppressed enough in two parameter regions [168, 169]. 
If axion-nucleon-nucleon interaction is weak enough, the axion cannot be effectively produced 
in the core of the supernova. Quantitatively, for fa ^ 10^ GeV, the axion flux can be small 
enough not to affect the cooling process [164-168]. On the contrary, if the axion interacts 
strongly enough, the mean free path of the axion becomes much shorter than the size of the 
core, and hence the axions cannot escape from the supernova. In this case, axion is trapped 
inside the so-called "axion sphere", and the axion emission is also suppressed. In this case, 
axions are emitted only from the surface of the axion sphere; this type of the axion emission is 
often called "axion burst". Quantitatively, for /q ^ 2 • 10^ GeV (or equivalently, nia ^ 3 eV), 
the axion luminosity from SN1987A is suppressed enough [164-169]. 

For /a ^ 2 • 10^ GeV suggested from the cooling of supernova, following the paper [169] 
we have another constraint from the detection of axions in water Cherenkov detectors. In 
this parameter region, the axion flux from the axion burst is quite sizable for its detection, even 
though it does not affect the cooling of SN1987A. If the axion-nucleon-nucleon coupling is strong 
enough, axions may excite the oxygen nuclei in the water Cherenkov detectors (^^0-|-a — )• ^^O ), 



^^In this context it is interesting to quote a remark from the paper by Hannestad S., Mirizzi A., Raffelt G.G., and 
Wong Y.Y.Y. [155]: "...In principle, fa ^ 10^ GeV is excluded by the supernova SN 1987A neutrino burst duration... 
However, the sparse data sample, our poor understanding of the nuclear medium in the supernova interior, and simple 
prudence suggest that one should not base far-reaching conclusions about the existence of axions in this parameter 
range on a single argument or experiment alone. Therefore, it remains important to tap other sources of information, 
especially if they are easily available" . 
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Figure 19: Astrophysical constraints on the axion mass nia from the cooling of the supernova, axion burst, cooling 
of the HB stars, the extragalactic background light [163] (squares), and the emission line in clusters of galaxies [162] 
(triangles). Shaded region is excluded. A purple star marks our result {ca-y-^ = 0.02, fa = 0.353 • 10^ GeV). The graph 
is inspired by [169]. 

followed by radiative decays of the excited state. If this process had happened, the Kamiokande 
detector should have observed the photons emitted from the decay of ^^O*. Due to the non- 
observation of this signal, /a ^ 3 • 10^ GeV is excluded [117]. 

Alternatively stated, the hadronic axion with the axion decay constant in the following range 
is still viable with all the astrophysical constraints^ [169] 



3 ^ ma ^ 20 [eV] 



(74) 



The mentioned consequences of the "axion trapping" effect may be illustrated as follows 
(Fig. 19) . Given that the value of the axion coupling to a photon in terms of the axion mechanism 
of Sun luminosity is ga-y ~ 7.07 • 10~^^ GeV~^, it is easy to write an equation for a straight line 
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106 GeV 

which marks a sort of a "border" (Fig. 19) between the allowed and forbidden values of the Ca-yy 
constant and the energy scale fa associated with the break-down of the U(l) PQ symmetry. 

Now everything is ready for discussion of the axionic contribution to the EBL. The axion 
with a mass nia = 17 eV has been termed "invisible" because it interacts very, very weakly. Its 
lifetime far exceeds the age of the Universe 



' a—f-'yy 



647r 



1439 • tun 



(76) 



where the age of the Universe tuniv ~ 4.34-10^'^ s [161, 170]. Notice that the lifetime of the axion 
is longer than the age of the Universe for nia = 17 eV and Ca-yy — 0.02 and hence primordial 
axions are still in the Universe. However, as we will see later, radiative decay of the axion may 
affect the background UV photons in spite of the long lifetime. 
The number (Na) of axions in a cluster (the mass M) is [171] 



^^Somewhat more conservative estimates performed by Turner [168] and Ressel [162] give the following limitations 
for the hadronic axion: 2 eV ^ nia ^ 5 eV and 3 eV ^ rua ^ 8 eV respectively. However, we are going to use the 
limitation (74) since the uncertainty of above-mentioned estimates is determined by the factor of ^ 3 [168]. 
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Msun \maj 

It seems to have gone unnoticed that this number may be so large that the cluster luminosity 

La = -^ (78) 

'a 

is easily measured [171]. In this connection we explore this possibility of observable photon 
luminosity caused by the cluster's axions decay. Let us remind that from now on we shall 
consider the diffuse extragalactic background (EBL) rather than the galaxy luminosity. 

According to the stated above, we shall start with examining the effect of dark matter in 
the form the light axions on the extragalactic background light. In dong so we follow the 
theoretical results by Overduin and Wesson [172], who assumed that the axions are clustered in 
Galactic halos with nonzero velocity dispersions and derived an expression for the intensity of 
the axionic contribution which describes the axion halos as a luminous element of a pressureless 
perfect fluid in the standard Friedman-Robertson- Walker universe. To go further and compare 
our predictions with observational data, we would like to calculate the intensity of axionic 
contributions to the EBL as a function of the wavelength Aq after the manner of [172]: 






Ao/(l+^)-A, 2 



dz 



0"A 



where Zf = 30 [172]; 



Xa = 24800A ( — ) (80) 



is a peak wavelength of the decay photons; 



2^A„ « 220A {'^\ 



ax = 2-Aa ~ 220A — (81) 

c \maj 

is the standard deviation of the Gaussian spectral energy distribution, for which the velocity 
dispersion Vc (for axions bound in galaxy clusters) rises to as much as 1300 km/s [162, 172]; 

f]. = 5.2- 10-3^2 (^) (82) 

is the present density parameter of the thermal axions; h is the Planck constant; Hq = IOO/iq km- 
s~^ ■ Mpc~^ is the Hubble constant, /iq = 0.75 it 0.15 [172] is the usual value of the Hubble 
constant expressed in units of 100 km • s^^ ■ Mpc^^; pcritfl = 1-88 • 10~^^ g ■ cm~^ is the present 
critical density; Vtmfi = ^a+^bar+^i> = 0.266±0.029 [173] is the present total density parameter 
of the axions (f^a), baryons (17;,^^ = 0.028 ± 0.012 [172]) and neutrinos {9.^ ^ 0.014 [174]); the 
expression for the decay lifetime of the axion decay into photon pairs was used in the following 
form: 

-«,,-^|-(i-5^-io'WC-(:i^)". C^li; (83) 

Evaluating Eq. (79) over 1500A ^ Aq ^ 20, OOOA with C = 0.03 and Zf = 30, we obtain 
the plots of I\{\q) shown in Fig. 20. Apparently, the theoretical fit of spectral intensity of the 
background radiation (83) produced by axion decays describes the known experimental data 
in the near ultraviolet and optical bands very well. They include data from several ground- 
based telescope observations (SS78 [175], D79 [176], BK86 [177]), sounding rockets (H77 [178], 
H78 [179]), Apollo-Soyuz mission (P77 [180]), the Pioneer 10 spacecraft (T83 [181]), the DIRBE 
instrument aboard the COBE satellite (H98 [182], WROO [183], COl [184]), S2/68 sky-survey 
telescope aboard TD-1 satellite (G80 [185]). 
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Figure 20: The spectral intensity /a(Ao) of the background radiation from decaying axions as a function of the observed 
wavelength Aq. The curves for the value of iria = 17 eV, C = Ca-yy/0-72 = 0.03 correspond to upper, median and lower 
limits on Hq. Also observational upper limits (solid symbols and a heavy line) and reported detections (empty symbols) 
over this waveband are shown. Experimental data depicted in red were not taken into account for the theoretical fit 
of the spectral intensity I\ (blue line). See explanations in the text. 

It is important to keep in mind that some of the experimental data were not taken into 
account at all during the construction of the theoretical fit of the spectral intensity Ix. For 
example, OAO-2 satellite (LW76 [186]), sounding rockets (J84 [187], T88 [188]), the Space 
Shuttle-borne Hopkins UVX experiment (M90 [189]), and combined Hubble Space Telescope - 
Las Campanas Telescope observations (B02 [190]), shown in red in Fig. 20. 

One of the reasons for not considering the data from LW76 [186], J84 [187], T88 [188], 
(M90 [189] is that these experiments involved measurements of certain parts of the sky only. 
This is a grave methodological disadvantage which, according to Gondhalekar [185], may lead 
to serious distortions of the true EBL value, because "...the individual observation were taken 
over different regions of the sky and cover different wavelength ranges. It should be noted 
that intensity of the observed inter-stellar radiation field shows significant variation, not only 
with galactic latitude but also with galactic longitude and these variations should be taken into 
account when comparing observation taken over different regions of the sky. Strictly, an accurate 
determination of the total interstellar radiation density requires integration of observations made 
over the whole sky". Experimental data from B02 [190] were dropped due to the other reason 
which is related to the fair criticism by Mattila [191], who casts doubt on the data calibration 
method and, consequently, on the accuracy of the obtained results. 

Let us now make a short comment regarding the quantity ("• For this purpose we write the 
effective Lagrangian density which describes the coupling ga-y of axions to photons: 



L. 



day 



a-yy 



F.^.F^'^c 



flU 



-Qa-yE ■ Ba. 



(84) 



where a is the axionic field, F is the electromagnetic field-strength tensor, F is its dual, and E 
and B are the electric and magnetic fields, respectively. The axion-photon coupling constant is 
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^'^^ - 2vr//'^^^ - 2^/J AT 3(1 + .) J' ^^'^ 

where Epq and A^, respectively, are the electromagnetic and color anomaly of the axial current 
associated with the axion field. Here z = mu/md is the u- and d-quarks masses ratio. 

Since the theoretical fit (Fig. 20) of the spectral intensity Ix of the background radiation 
from decaying axions was performed for the value C = 0.03, the constant Ca-y-y, according to (83), 
will be 



L^a-y-y 



0.72C 



EpQ 2(4 + z) 



N 3(1 + z) 



0.02. (86) 



This value is extremely important, because according to [192], it is a strong indicator of (a) 
the effects of axion emission on the evolution of helium burning low-mass stars [193], (b) the 
effect of decaying relic axions on the diffuse extragalactic background radiation [162, 163]. That 
is the effects (a) and (b) provide limit to the axion-photon coupling [192]. It can be shown, that 
a combined action of effects (a) and (b) with z = niu/md = 0.56 and fa — 10^ GeV leads to 
the following important limitation: 



Ca77 ^ 77^7^^7777 = 0-02. (87) 



A 

10^ GeV 

The exact match of Eqs. (86) and (87) reflects a remarkable fact that the properties of the 
axion studied {rua = 17 eV, fa = 0.353- 10^ GeV) conform to conditions (a) and (b) completely. 
Moreover, they satisfy the conditions (see [192] and Fig. 19) of the axion's effect on the neutrino 
burst from SN 1987A (the so-called "trapping regime" in which the axion emission would not 
have a significant effect on the neutrino burst [164-168]) and the effect of axions emitted from 
SN1987A on the Kamiokande II detector, coming from a condition of absence of a large signal 
at the Kamiokande detector [117]). 

Noteworthy, the axion under question {ma = 17 eV, Cayy = 0.02) may be put into the 
thermally- produced hadronic axions class (see Fig. 17 and Fig. 18), and it may also be considered 
as a real candidate for dark matter at the same time. Here is why. 

It is usually assumed that dark matter in standard cosmology models is produced during 
the radiation-dominated era. If the axion mass rua ^ 10~^ eV, axions were produced thermally, 
with cosmological abundance 

^ahl = -^ (^ , (88) 

130 eV \g*s,Fj 

where g*s,F is the effective number of relativistic degrees of freedom when axions freeze out 
of equilibrium. Turner showed [164, 165] that the vast majority of these would have arisen in 
the early Universe via thermal mechanisms such as Primakoff scattering and photo-production. 
The Boltzmann equation can be solved to give their present comoving number density as na = 
{830/ g^:S,F) cm~^ [162], where g^:S,F ~ 15 counts the number of relativistic degrees of freedom 
left in the plasma at the time when axions "froze out" of equilibrium. The present density 
parameter 0^ = narua/ Pcritfi of thermal axions thus leads to the expression (82) whence it 
follows 

0.11 ^ Sla ^ 0.25. (89) 

Here we have taken 0.6 ^ ho ^ 0.9 as usual. This is comparable to the density of dark 
matter. 

At the same time, there are axion constraints in the so-called nonstandard thermal histories. 
As it was shown in [194] , the most intriguing one among the well-known non-standard thermal 
histories {low-temperature reheating and kination cosmologies) is the LTR cosmology which 
allows for the fact that there is currently no direct evidence for radiation domination prior to 
big-bang nucleosynthesis [195]. According to this scenario, radiation domination begins as late 
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Figure 21: (a) Upper limits to the hadronic axion mass from cosmology, allowing the possibility of a low-temperature- 
reheating scenario. The green region shows the region excluded by the constraint flah^ < 0.135 as a function of 
the reheating temperature Trh- The red region shows the additional part of axion parameter space excluded by 
WMAPl/SDSS data. At low reheating temperatures, upper limits to the axion mass are loosened. For Trh ^ 170 MeV, 
the usual constraints are recovered. Adapted from [194]. 

(b) Estimated improvement in the accessible axion parameter space from including more precise measurements of 
the matter power spectrum (the region bounded by the yellow line), corresponding to LSST [196, 197], or from mea- 
surements of clustering on smaller length scales, corresponding to Lyman-a forest measurements (the region bounded 
by the blue line) [198]. The red region indicates the parameter space excluded by WMAPl/SDSS measurements. 
Adapted from [194]. Our result (a red star) has the coordinates {Trh = 30 MeV; rua = 17 eV}. 



as 1 MeV, and is preceded by significant entropy generation. Thermal axion relic abundances 
are then suppressed, and cosmological limits to axions are loosened. However, for reheating 
temperatures Trh ^ 35 MeV , the large-scale structure limit to the axion mass is lifted (see 
Fig. 21). The remaining constraint from the total density of matter is significantly relaxed in 
this case. Constraints are also relaxed for higher reheating temperatures. It turns out that 
axions will be produced thermally, with nonstandard cosmological abundance 



nth 1,2 



ViT^k 



nia 







10 



liTrh/Tp) 



(90) 



130 eV \g*s,F, 

where Tp is the decoupling temperature of axions. 

Using the abundance Q"*'^ normalized by its standard value 17^ as a function of the reheating 
temperature, obtained by Grin et al. (Fig. 4 in [194]) it is easy to derive for the mass nia = 17 eV 
at Trh = 30-35 MeV (see Fig. 21) 



^nth 



-f {Trh /Tp) ^ 0.2 for Trh = 30 MeV. 



Hence, 



0.022 ^ nf' < 0.05. 



(91) 



(92) 



which is rather small and comparable to the baryonic density value 0.027 ^ 0^ ^ 0.040 [161] as 
opposed to the standard case of (89). 

Which one of them is preferable when it comes to description of real physics - the standard 
expression (89) for the axion dark matter density or the non-standard one (92)? Although 
there is no answer to this question nowadays, one may still hope that, according to [194], future 
limits to axions in the standard radiation-dominated and LTR thermal histories may follow from 
constraints to their contribution to the energy density of relativistic particles at T ~1 MeV. 
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This is due to the fact that a comparison between the abundance of ^He and the predicted 
abundance from standard big-bang nucleosynthesis (SBBN) places constraints to the radiative 
content of the Universe at T ~1 MeV [199]. As the paper [194] notes, this can be stated as a 
constraint to the effective neutrino number {Ny^), because at early times, axions contribute 
to the total relativistic energy density (through N^ ), and thus constraints to ^He abundances 
can be turned into constraints on nia and Tj.^. 

It is known that in terms of the baryon-number density n^, the primordial He abundance is 
characterized by the expression Yp = 4n/fe/'i-6- In order to translate measurements of primordial 
^He abundance Yp to constraints on nia and Tj-h we use the scaling relation by Grin et al. [194] 

4S 

AN'Jf{ma,Trh) = N'Jf - 3 = y {(6-25 • AYp + !)-!}, (93) 

which describes the relation between the deviation AN^^ for the given values of nia and Trh 
and the deviation AYp of the primordial '^He abundance. Here the deviation AYp is thought 
of as a deviation from the value Yp = 0.2487 it 0.0006 accepted for SBBN-predicted primordial 
abundances [200]. 

According to the calculations in [194], for the values rua = 17 eV, Trh = 30 MeV and 
the effective neutrino number A'^^ = 3.44 (Fig. 7 in [194]), the deviation ANf, is 0.44. 
Substitution of this value into Eq. (93) gives the following value for the AYp of the primordial 
'^He abundance: 

AYp = 0.0044, (94) 

Let us present some data for the sake of comparison. One careful study gives the value 
Yp = 0.2565 ± 0.0010 {stat) ± 0.0050 [syst] from 93 Hn regions [201], leading to the sensitivity 
limit Ny^ = 3.61. Interestingly enough, the deviation ANt = 0.61 may be a sign of the 
higher mass axions existence, since according to [192], for sufficiently high masses, the axionic 
contribution saturates to AN^ = 4/7 at high reheating temperatures. At the same time 
the Planck satellite is expected to reach AYp = 0.013 [202], yielding sensitivity of Nu — 
4.04, while CMBPol (a proposed future CMB polarization experiment) is expected to approach 
AYp = 0.0039, leading to the sensitivity hmit K^-'' = 3.30 [194]. 

It makes it clear that the answer to a question about the preference of one expression 
(standard (89) or non-standard (92)) for the axion dark matter density over another may be 
searched for only in direction of a sharp increase in observations precision. According to [203] , 
measuring CMB temperature and polarization with cosmic variance accuracy would allow to 
constrain Yp to within 1.5%, or AYp ~ 0.0036 (assuming flatness). Such an ideal measurement 
would be able to discriminate between the BBN-guided, deuterium based helium value and the 
current lowest direct helium observations. In other words, "...if the CMB-determined helium 
mass fraction turns out to be as high as suggested by SBBN calculations combined with the 
observed deuterium abundance, this could indicate a systematic error in the present direct 
astrophysical helium observations. Alternatively, if the CMB could independently determine 
the helium value with sufficient precision to confirm the present low helium value coming from 
direct observations, then this would be a smoking gun for new physics" [203]. For example, one 
could imagine sterile neutrinos appearing within the nonstandard BBN scenarios, which would 
agree with present observations of r/io = W^^{nb/n^), while having a low helium mass fraction. 
To put it differently, in our opinion, in spite of the future possible constraints to axions and 
low-temperature reheating from the helium abundance and next-generation large-scale-structure 
surveys, it is the appearance of the sterile neutrinos that may effectively solve the problem of 
missing dark matter in the framework of LTR cosmology. 

And finally, returning to our 17 eV axion "caught" in the extragalactic background, we 
may say that regardless of the scenario which provides a significant portion of the dark matter, 
experimentally observed invisible axions (Fig. 20) must have rest masses in the " semi- visible" 
range (1500A - 20000A) where they do contribute significantly to the light of the night sky. In 



41 



this sense one may think of our findings as a result of the axion mechanism of Sun luminosity 
and the "...nature's most versatile dark-matter detector: the light of the night sky" [172]. 

6 Relic axion-like archion and cosmic infrared 
background 

It is well known that a direct measurement of the EBL and, particularly, the cosmic infrared 
background (CIB) consists in observing the cumulative emission from various pregalactic ob- 
jects, protogalaxies, galaxies, cosmic explosions and decaying elementary particles (including 
dark matter particles) throughout the evolution of the Universe and therefore one can provide 
important constraints on the integrated cosmological history of star formation [204, 205] and 
control the second most important contribution to the cosmic electromagnetic background after 
the Cosmic Microwave Background generated at the time of recombination at a redshift around 
1000 [205]. This background is expected to be composed of three main components [205]: 

• the stellar radiation in galaxies concentrated in the ultraviolet and visible with a redshifted 
component in the near InfraRed (IR); 

• a fraction of the stellar radiation absorbed by dust either in the galaxies or in the inter- 
galactic medium; 

• the radiation from active galactic nuclei (a fraction of which is also absorbed by dust and 
reradiated in the far-IR). 

Its detection is a subject of great scientific interest and the main purpose of the Diffuse 
Infrared Background Experiment (DIRBE) on the Cosmic Background Explorer (COBE) space- 
craft [182-184]. These studies resulted in upper limits on the EBL in the 1.25-100 fim region, 
and in the detection of a positive isotropic signal at 140 and 240 fim. 

However, there are serious problems with interpretation of some DIRBE data at far-IR wave- 
lengths. For example [204], the energy sources could either be yet undetected dust-enshrouded 
galaxies, or extremely dusty star-forming regions in observed galaxies, and they may be respon- 
sible for the observed iron enrichment in the intracluster medium. Although there is currently no 
compelling need to invoke non-nuclear energy sources to explain the COBE data, their potential 
contribution to the observed EBL cannot be ruled out. It leads to a conclusion that the exact 
star formation history or scenarios required to produce the EBL at far-IR wavelengths cannot 
be unambiguously resolved by the COBE observations and must await future observations [204]. 

The other type of problems is revealed by the interpretation of the DIRBE data at near-IR 
wavelengths. For example, in their studies of EBL at near-IR wavelengths Cambresy et al. [184] 
obtain a significantly higher cosmic background than integrated galaxy counts (3.6ib0.8 kJy-sr^^ 
and 5.3±1.2 kJy-sr~^ for 1.25 jim and 2.2 //?7i, respectively), suggesting either an increase of the 
galaxy luminosity function for magnitudes fainter than 30 or the existence of another contribu- 
tion to the cosmic background from primeval stars, black holes, or relic particle decay. However, 
models predict other possible contributions to the background at these wavelengths [206] such 
as a burst of star formation either in primeval galaxies or in Population HI stars {z ~ 10), very 
massive black holes (accreting from a uniform pre-galactic medium at z ~ 40), massive decaying 
big bang relic particles {z ~ 300). In this regard, Cambresy et al. [184] make a natural conclu- 
sion that the new constraints in the near-infrared should encourage revisiting the importance 
of those contributions to the CIB in cosmological models. 

In this connection we consider the effect of dark matter in the form of 2 eV relic particle on 
the extragalactic background light at near-IR region (Fig. 22). It was obtained on the basis of the 
intensity (79) of axionic contributions to the EBL calculation as a function of the wavelength 
Aq. The theoretical fit (green curves in Fig. 22) of the spectral intensity of the background 
radiation (79) produced by axion decays, apparently, describes the known experimental data for 
the near-infrared band very well. They include data from ground-based telescope observations 
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Figure 22: The spectral intensity /a(Ao) of the background radiation from decaying axions as a function of the observed 
wavelength Aq. The curves for values rua = 17.0 eV, ( = Caj^/0.72 = 0.03 and TOq = 2.0 eV, C, = 9.0 correspond 
to upper, median and lower limits on /ig. Also observational upper limits (solid symbols) over these wavebands are 
shown. 



(BK86 [177]) and the DIRBE instrument aboard the COBE satellite (H98 [182], WROO [183], 
COl [184]). 

Let us consider some properties of such relic particle with the 2 eV mass. It may be imme- 
diately stated that such particle cannot play the role of a sterile neutrino with the mass of the 
order of a few eV [207]. And the reason is the following. 

It is known that if neutrinos are massive and if the mass eigenstates are not degenerate, 
then it is possible to have a radiative decay of the form i/^ — )■ i/ + 7. According to [208], this 
gives for the decay rate of Majorana neutrinos 






^1/7 



= 5.52 • 10"^^ 



^sin^g^ 



\kev) "" ' 



(95) 



where ms is the mass eigenstate most closely associated with the sterile neutrino, and 9 is 
the mixing angle between the sterile and active neutrino. The decay of a nonrelativistic sterile 
neutrino into two (nearly) massless particles produces a line at the energy E'j = ms/2. 

Obviously, in case of the decay reaction i/^ — )• z^ + 7, Eq. (79) can be generalized to another 
relic Vs by multiplying it by 



/a(Ao 



1 Q r* , 

-L i^u a— >-77 

9 O T* ' 



(96) 



where the 1/2 is a number of photons produced in each i/g decay, Vt^, ^ 0.014 is the present total 
density parameter of the neutrinos. 



>!/7 



(r: 



us^iyy/ 



(1.7.10-t.™.)C-fS) 



\keVJ 



(97) 



where 
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Using the generalization (96), it is rather easy to show that a theoretical fit of spectral inten- 
sity I\ of the background radiation is several orders of magnitude lower than the experimental 
data in near-infrared band under any reasonable values of the mixing angle sin^ 29. 

This relic particle with the 2 eV mass may be supposed to belong to a class of axion-like 
particles with rather exotic properties. One of them may originate from the fact that in the 
distant past their birth in the electromagnetic field, or in other words interaction with a photon, 
was not suppressed for some reason, while at the present time there are simply no conditions 
suitable for their birth. To put it differently, the mentioned conditions of such particles' birth 
must be completely suppressed nowadays in order our axion mechanism of Sun luminosity to 
be possible. 

In our opinion, such axion-like particle, provided that it exists, has the properties most 
similar to those of the so-called archion [209-215]. An archion may be very similar to a hadronic 
axion with highly suppressed interaction with leptons under certain conditions. Let us discuss 
briefly some of its properties below. 

As is generally known, in all the models of an invisible axion this particle appears as a 
Goldstone boson connected with the phase of a complex SU(2)xU(l) singlet Higgs field. The 
axion coupling to the gauge bosons appears in these models after the U(1)pq symmetry violation 
by means of a mechanism, specified by the non-vanishing color anomaly U(1)pq - SU(3)c - 
SU(3), [216]. 

In its most generic form the Lagrangian of axion interaction with fermions (quarks and 
leptons) and photons is 

L = Ca/sa ■ FQ,(sin 0ai3 + ij5 cos 6'q,^)-F/3 + ga^aF^t.F^" , (99) 

where a, /3 = 1, 2, 3 indices denote the generation of fermions F, and the constants OajS, 

Caf3 OC /-I (100) 

and 

9ay CC /-I (101) 

depend on the axion model chosen. 

A model of an archion [209-215] arose from a model of horizontal unification, the basis of 
which is expounded in a monograph by Khlopov [216]. This theory includes the global 1]{1)h 
symmetry, the spontaneous breaking of which leads to prediction of a Goldstone boson of the 
invisible axion type. Such boson called "archion" by authors of [211-214] has the flavor non- 
diagonal as well as the flavor diagonal coupling to fermions. 

The global \]{1)h symmetry in horizontal uniflcation may be identified with the Peccei-Quinn 
symmetry U(1)pq [28-35], which is due to the fact of triangle anomaly existence in axial currents 
\J{1)h interaction with gluons [216]. 

In the simplest variant of the horizontal unification (the gauge symmetry 

S{U)c <E) SU{2) U{1) <E) SU{3)h U{1)h (102) 

with a minimal set of heavy fermions) the anomaly is compensated, and the archion remains 
almost massless. The interaction of the archion with photons is absent because of the parallel 
compensation which corresponds to the current-photon-photon anomaly. 

On the other hand, according to [216], within any realistic extension of the horizontal uni- 
flcation model up to the grand uniflcation symmetry, for example, within the extension to the 
SU(5)/f(8)SU(3)// symmetry, there is no compensation because of the extra heavy fermions, so 
that the archion appears to be similar to a hadronic axion with strongly suppressed interaction 
with leptons. 
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In the framework of the archion model, the archion-photon couphng constant ga-y which 
appears in Eqs. (99) and (101) has the following form [216]: 



Say 



a a AcZ^I"^ 



27r fa"'' 47r/,(l + z) 



Aem 2(4 + Z) 

Ac 3(1 + z) 



(103) 



Here as usual z = niu/md is the u- and d-quarks mass ratio, fa is the energy scale associated 
with the breakdown of the U(l)// (or equivalently, the U(1)pq symmetry), Ac and Aem are the 
color and electromagnetic anomalies respectively. 

Taking into account the value of the Ca^y = 0.72(^ constant (see (83) and Fig. 22), which is 
~6.5 for the archion mass of 2 eV, it is easy to estimate the archion-photon coupling constant 
from (103) and (10): 

5*^-2.5-10-9 GeV-^. (104) 

Obviously, if the archion indeed exists and has such a high archion-photon coupling constant, 
it must be of the relic origin only in a sense that the conditions of such particles birth must 
be completely suppressed nowadays. Otherwise, our axion mechanism of Sun luminosity would 
not be possible, like we pointed out earlier^®. The value of this constant is also very important 
because according to [192], it is a strong indicator of (a) the effects of axion emission on the 
evolution of helium burning low-mass stars [193] , (b) the effect of decaying relic axions on the 
diffuse extragalactic background radiation [162, 172]. In other words, an archion with the 2 eV 
mass, characterized by the archion-photon coupling constant (104), must be relic in order not 
to violate the known limitation (87) which is a consequence of combined action of effects (a) 
and (b). 

For all invisible axion models, including the archion model, Lagrangian of its interaction with 
nucleons has the same form [215]. Given (104), it is not hard to estimate the archion- nucleon 
coupling constant for the 2 eV mass archion (Fig. 22) using (29)-(31) and (10): 

51 ^ 5.6 • 10-8. (105) 

If the archion exists with such low value of the archion- nucleon coupling constant, it obviously 
must have only relic origin in order not to violate the SN1987A limit (3-10~'^ ^ gan ^ 10^^ [116, 
117]). 

This rises the question of why would Nature need a relic axion-like archion in addition to 
an ordinary hadronic axion. Curiously enough, the answer suggests itself and is related to the 
evolutionary formation of the visible large-scale cosmological structure against the background 
of the invisible "dark" structure. 

Apart from the details of this scenario, our axion-like particles may be among the primary 
participants of this creative action, if they do exist. Let us note however that in spite of the fact 
that the cosmological structure formation is provided by the single kind of hidden-mass particles 
- axions, the "axion liquid" turns out to be two-component. The visible structure in the form 
of galaxies and superclusters is formed by the shortwave component of the 17 eV axions density 
perturbations spectrum, while the relic thermal 2 eV archions background plays an important 
role in sub-shortwave density perturbations spectrum evolution. The latter includes, in partic- 
ular, the massive halos formation beyond the visible parts of the galaxies. Following [216], it is 
worth mentioning that the so-called phase-space argument by Tramaine and Gunn [217], which 
gives rise to a known limit on particles mass in the halo, may be substantially weakened or even 
omitted for the Bose gas [218-221]. 

If all said above is true, it is more than sufficient for a "sensible" existence of axions, and 
particularly, a relic axion-like archion. On the other hand, except for the indirect indication of 
axions existence in the form of Fig. 22, we have no other - direct - arguments. In order to fill 



^*It is necessary to note that regardless of whether one assumes the existence of the axion mechanism of Sun 
luminosity or not, such high value of the archion-photon coupling constant (104) is forbidden by the DAMA experiment 
observations [151] presented in Fig. 17a and 17b. 
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this gap a little, let us try to substantiate the quantitative relations between the experimental 
spectral intensities of the background radiation from decaying axions and archions shown in 
Fig. 22 theoretically on the basis of a simplified cosmological model which takes into account 
the processes of axion dark matter radiative decay. 

6.1 Decaying axion and relic archion as two components of 
luminous dark matter 

Let us consider the flat (for simplicity) Universe after the recombination consisting of the follow- 
ing components: relic radiation, usual (light, visible) matter, dark matter (presented by axions 
and archions), non-relic radiation (resulting from their decay) and dark energy (described by 
the cosmological constant A), then for the corresponding metrics 

ds'^ = c^dt^ - a^{t){dx^ + dy"^ + dz^), (106) 

where a{t) is the scale factor, the first Friedmann equation reads 

^^ = AcToo + A, (107) 

where a dot denotes the derivative with respect to t, k = SttGn/c^ {Gjy is the Newtonian 
gravitational constant) and the 00 covariant component of the total energy-momentum tensor 
reads 

Too = ^rr + ^i/m + £a2 + £al7 + £r2 + ^rlT, (108) 

where, in their turn, err,um,ai,ai7,r2,rn denote energy densities of relic radiation (assumed to be 
independent of all other components), visible matter (also assumed to be independent), archions 
(with the mass ma2 = 2 eV and the disintegration constant Xa2 = 3.88 • 10~^^ s, assumed to be 
completely nonrelativistic), axions (with the mass man = 17 ^.V and the disintegration constant 
Aai7 = 1.91 • lO^^^s^^, also assumed to be completely nonrelativistic), radiation resulting from 
archion decay and radiation resulting from axion decay respectively. 

Here some additional comments should be made. First, the disintegration constants are 
estimated simply as A = 1/r* on basis of the formula (83). Second, we assume both dark matter 
components (archions and axions) completely nonrelativistic, in other words, we completely 
neglect their velocities (or temperatures). This assumption is simultaneously in agreement and 
disagreement with [172], where the formula (185) (which is equivalent to (79) in our text) 
contains simultaneously the scent of the standard cosmological ACDM-model without matter 
velocities in the denominator of the integrand and the scent of the velocity dispersion in its 
numerator. It seems that this fact does not mean that there is a self-contradiction in the 
formula (185), because the numerator may be much more sensitive to the local velocity of dark 
matter than the denominator to the global one. However, in the subsequent analysis we shall be 
interested in the evolution of the average non-relic radiation energy density without taking into 
account its frequency distribution and the nonzero dark matter temperature (for simplicity). 

Substituting (108) into (107) and introducing the standard portions 

2 2 

^ = ^^m, ^^ = ^^' (109) 

where the subscript (0) corresponds to the current moment of time t = (this value may be 
chosen without loss of generality) and H = d/a is the Hubble parameter, we obtain 

~ 1 ~ -^0 ( ^rr i" ^um 1" ^a2 h 

Oy V err(O) ^um{0) ^a2{0) 

+nal7^^"^+^r2^^+^rl7^^^^+^A). (HO) 

eal7(0) ^r2(0) ^rn{0) J 
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Here Err ~ 1/a^, ^um ~ l/o^^, consequently, as usual, 
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where qq is the current value of the scale factor: a(0) = ao- Further, 
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Finally, combining first and second Friedmann equations, we come to the following equations: 



1 1 

d [{ea2 + er2)a^] + -er2d{a^) = 0, d [{san + ern)a^] + -ernd{a^) = 0. 



From the first one we immediately get 
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Similarly, 



Eriy = — I eal7(0)aoAal7 



/ a{t) exp{—Xai7t)dt -\- e. 



rl7(0)Oo 



(118) 



The substitution of (111), (113), (117) and (118) into (110) gives the equation defining 
the function a{t) satisfying the condition a{t) = oq. Taking into account that the values of 
both disintegration constants Aa2 and Aai7 are extremely small and the radiation energy density 
decreases faster than the nonrelativistic matter energy density when a increases, one can neglect 
all (i.e. both relic and non-relic) radiation contributions and get 
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Introducing the dimensionless quantities 



T TT . ~ "^ \ Aa2 T ^al7 
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from (112) we obtain 
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Figure 23: The numerical solution of Eq. (121) (the blue firm line) is indistinguishable from its numerical solution for 
Aa2 = and Aqiy = (the yellow dashed line). 
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According to the recent observations [222], we consider the following values: ^um = 0.046, 
^a2 + ^ai7 = 0.236, f^A = 1 — ^um — (^a2 + ^au) = 0.718. For Qa2 and fiair we use the values 
0.016 and 0.220 respectively. Again, due to extreme smallness of Xa2 and Xan (in the presence 
of nonzero A) they do not affect noticeably the dependence a{t). The following graphs confirm 
this statement: 

From (117) and (118) we obtain respectively 



id{i) exp{—\a2i)di + 
td{t) exp{—Xani)di + 
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Let us introduce the convenient functions 
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(124) 



(125) 



For a = we get fa2 = —9.26 • 10~^ and fan = —4.56 • 10~^. These negative values may be 
completely compensated by the proper choice of £,-2(0) and £^17(0) ^^ order to give £,-2 = and 
£ri7 = at the same moment of time (when a = 0). This has clear physical sense: we assume 
that at the recombination moment the non-relic radiation is absent (all radiation, present at that 
moment of time, may be considered as the relic one). From this requirement we immediately 
obtain £,,2(0) = 9.26 • 10-^£„2(o) and £,.17(0) = 4.56 • 10-^£ai7(o). 

Now let us also depict the helpful graphs of the ratios £r2/£a2 and £ri7/£ai7 as functions of 
t (Fig. 24a, red and green lines respectively): 



£^2 exp(Aa2t) I I /'wrx f T rx ,7 , ^^-2(0) 
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Figure 24: Graphs of the ratios SriJEa^ (red), eri7/£ai7 (green) (a) and £r2/£ri7 (b) as functions of time t. 
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£al7(0) 

In particular, for t = we get er2l^a2 = 9.26 • 10~^ and Srvj jsavj = 4.56 • 10^^, as it should 



be (Fig. 24a). 

Let us also depict the helpful graph of the ratio Er^l^rii as a function of t (Fig. 24b). 
The theoretical estimate of the ratio Eril^rYi is, apparently, equal to 



£r2 
ErYI 



0.15. 



(128) 



theory 



On the other hand, it is rather easy to estimate the experimental er2/£ri7 ratio using the 
data shown in Fig. 22 and the expression for the spectral intensity I\{\q) of the background 
radiation from decaying axions: 
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(129) 



Curiously enough, the theoretical estimate (128) and the experimental one (129) agree fairly 
well, which increases one's optimism as to the essence of the results obtained in the current 
section. At the same time, our estimates require additional thorough verification, especially 
when it comes to the experimental justification of the validity and reliability of the data on 
the extragalactic background light in near- ultraviolet, optical and near- infrared bands (1500- 
lOOOOA). Indeed, it is a very important step to perform, because otherwise our unusual results 
may seem more like just the luckily guessed rules of calculation, which do not reflect the actual 
nature of things. 

7 Summary and Conclusion 

In the present paper we present a self-consistent model of the axion mechanism of Sun luminosity 
and solar dynamo ~ geodynamo connection, in the framework of which we estimate the values 
of the axion mass {nia ~ 17 eV) and the axion coupling constants to photons (g'a^ ~ 7.07 • 



-11 



10 ^^ GeV ), nucleons (^an ~ 3.20- 10 ) and electrons ((^ae ~ 5.28 • 10 ). Their verification 
on the basis of the model results comparison with the known experiments is also provided, 
including the CAST-, CUORE- and XMASS-experiments. 
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In order to explain the solar-terrestrial magnetic connection we propose the axion mechanism 
explaining the Sun luminosity physics and "solar dynamo - geodynamo" connection, where the 
total energy of axions, which appear in the Sun core, is initially modulated by the magnetic 
field of the solar tachocline zone due to the inverse coherent Primakoff effect and after that 
it makes its way towards the Earth where its "iron" component containing the 14.4 keV solar 
axions is resonantly absorbed inside the nickel-iron core of the Earth. It results in the fact that 
the variations of the axion intensity play a role of an energy source and a modulator of the 
Earth magnetic field. In other words, the solar axion mechanism is not only responsible for 
formation of a thermal energy source in the liquid core of the Earth necessary for generation 
and maintenance of the Earth magnetic field, but unlike other alternative mechanisms [15], 
naturally explains the cause of the experimentally observed strong negative correlation of the 
magnetic field in the tachocline zone of the Sun and the magnetic field of the Earth. 

It is necessary to note that obtained estimations can't be excluded by the existing experi- 
mental data (see Fig. 17a and Fig. 18a) because the effect of solar axion intensity modulation 
by temporal variations of the toroidal magnetic field of the solar tachocline zone discussed 
above was not taken into account in these observations (see Fig. 17b and Fig. 18b). On the 
other hand, the obtained estimates for the axion-photon coupling and the axion-nucleon cou- 
pling cannot be ruled out by the existing theoretical limitations known as the globular cluster 
star limit {ga-y < 6 • 10~^^ GeV~^) and the red giant star limit {gae < 3 • 10~^^) [105], since 
these values are highly model-dependent. It actually means that the axion parameters ob- 
tained in the present paper do not contradict any of the known experimental and theoretical 
model-independent limitations. 

Let us give some major ideas of the paper which formed a basis for the statement of the 
problem justification, and the corresponding experimental data which comport with these ideas 
either explicitly or implicitly. 

7.1 Axion mechanism of Sun luminosity 

One of the key ideas behind this mechanism is the effect of 7-quanta channeling along the 
magnetic flux tubes (waveguides) in the Sun convective zone (Fig. 8), which may be represented 
by the 7-quanta channeling in the periodical structure (Fig. 6) in the particular case. A low 
refraction (i.e. a high transparency) of the thin magnetic flux tubes is achieved due to the 
ultrahigh magnetic pressure (see (16)), induced by the magnetic field of about 200-400 T (Fig. 9a, 
adapted from [79]). It is noteworthy that although such strong magnetic fields have never been 
used for the explanation or interpretation of the simulation results such as the stability analysis 
of tachocline latitudinal differential rotation and coexisting toroidal band using an MHD analog 
of the shallow- water model [79], they have always been present implicitly in the very same 
simulation results, as our analysis shows. Fig. 9a provides an illustrative example where a 
hidden part of the "latitude - magnetic field of the overshoot tachocline zone" dependence 
(which is absent on the original plot in Fig. 11 of [79]) is added. A right to exist and physical 
validity of this part is confirmed by the bare fact of its direct correspondence (Fig. 9a, b) with the 
real X-ray image of the Sun in its active phase (Fig. 9b) obtained in the experiments performed 
with the Japanese X-ray telescope Yohkoh (1991-2001) (adapted from [14]). 

The direct experiments on monitoring the space-time evolution of the magnetic flux tubes, 
rising from the deep layers of the convective zone of the Sun (Fig. 8b, adapted from [67]) proved 
the existence of the "hollow" ideal magnetic waveguides for 7-quanta. These tubes cross the 
photosphere and form the solar active regions such as sunspots which are the sources of X-rays 
(Fig. 8d, adapted from [71]). There are also some theoretical results substantiating the effect 
of anchoring magnetic flux tubes in the tachocline (Fig. 8c, adapted from [70]). 

Finally, the most impressive evidence of the axion mechanism of Sun luminosity are the solar 
images (Fig. 13, adapted from [14]) taken at photon energies from 250 eV up to a few keV from 
the Japanese X-ray telescope Yohkoh (1991-2001), which depict the solar X-ray activity during 
the last maximum of the 11-year solar cycle (Fig. 13b). It is hard to imagine another model or 
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considerations which would explain such anomalous X-ray radiation distribution over the active 
Sun surface just as well. One should also keep in mind that the alternative mechanism of axion 
transformation into 7-quanta in the magnetic field of sunspots or other solar active regions is 
completely ruled out, since there are no signs of an X-ray bright spot at the disk center (see 
Fig. 13a and b), which should otherwise be observed according to [223]. 

7.2 Invisible axions and Solar Equator effect 

According to the axion mechanism of Sun luminosity, a part of axions that pass through the 
tachocline zone near the equator and poles (Fig. 6 and Fig. 13) is not converted into 7-quanta 
by the Primakoff effect because of the magnetic field vector collinearity to the axion momentum. 
It is the equatorial part (Fig. 13) of invisible axions that reaches the Earth, where its "iron" 
component (i.e. the 14.4 keV solar axions) is resonantly absorbed in the Earth core. The 
energy of the solar axions supplied to the Earth core in this way plays a role of a trigger for 
the generation and maintenance of the geomagnetic field, thus giving birth to the effect of the 
steady anticorrelation between the variations of the Solar magnetic field and geomagnetic field 
(Fig. 3). 

In this sense, it is appropriate to make a short remark related to the anticorrelation between 
the solar and terrestrial magnetic fields variations. The small variations of TSI (ATSI) are 
apparently produced by the solar magnetic field variations. At the same time, these variations of 
the solar magnetic field drive the equatorial sector width variations (Fig. 13), and consequently, 
the "equatorial" axion flux (see Fig. 12) 

A^^„ = 0.05TAr5//(r5/). (130) 

In other words, the "equator" effect is not only the source of "invisible" axions, it also 
modulates their intensity inversely proportional to the solar magnetic field change (see (130)), 
thus maintaining the inverse correlation between the solar magnetic field and the "invisible" 
axions flux (Fig. 12). The latter is also the cause of the inverse correlation between the solar 
and terrestrial magnetic fields variations (Fig. 3). 

7.3 Axion mechanism of the solar dynamo — geodynamo 
connection 

This mechanism is responsible for the physical (axionic) nature of the steady anticorrelation 
between the variations of the solar magnetic field and the geomagnetic field (the Y-component) 
which is directly proportional to the westward drift of magnetic features (Fig. 4, Fig. 5) at the 
measurement points (Western Europe and Australia). Figuratively speaking, in this case the 
Y-component of the Earth magnetic field acts as a "measuring instrument" which tracks the 
influence of the "equatorial" 14.4 keV solar axions on thermal and magnetic processes in the 
liquid core of the Earth. The physical scenario for this may be following. 

The "iron" solar axions that are resonantly absorbed in the Earth core, activate the vertical 
background motion along the gravity force. This motion, in its turn, "pushes" the temperature 
gradient to the bottom of the liquid core more or less heavily (depending on the total energy of 
axions), thereby changing the temperature profile in the convective medium of the Earth core. 
The change in the temperature profile also changes the thermal conditions near the nuclear 
georeactor situated at the boundary of the liquid and solid core of the Earth. Since the power 
output of such reactor is proportional to temperature in the range of 3000-50000° typical for 
the Earth liquid core, it means that the variations of the Earth core temperature generated 
by the mechanism of solar dynamo-geodynamo connection induce the corresponding variations 
of the nuclear georeactor thermal power. If the georeactor hypothesis is true, the fluctuations 
of georeactor thermal power, induced by the variations of the absorbed axions energy in the 
Earth's core, can influence the Earth global climate in the form of anomalous temperature jumps 
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in the following way. Strong fluctuations of the georeactor thermal power can lead to partial 
blocking of the convection in the liquid core [118-120] and the change of the angular velocity of 
liquid geosphere rotation, thereby changing the angular velocities of the Earth mantle and crust 
by virtue of the total angular momentum conservation law. It means that the heat, or more 
precisely the dissipation energy caused by friction of the Earth surface and lower atmosphere, 
can make a considerable contribution to the total energy balance of the atmosphere and thereby 
influence on the Earth global climate evolution significantly [118-120]. 

On the other hand, it is clear that understanding of the mechanism of the solar-terrestrial 
magnetic correlation can become the clue of the so-called problem of solar power pacemaker 
related to possible existence of some hidden but crucial mechanism of the Sun's energy influence 
on the fundamental geophysical processes. It is interesting that the "tracks" of this mechanism 
have been observed for a long time and manifest themselves in different problems of solar- 
terrestrial physics, in particular, in climatology where the mechanisms, by which small changes 
in the Sun's energy output during the solar cycle can cause change in the weather and climate, 
have been a puzzle and the subject of intense research in recent decades. 

If we also add the fact that the dissipation energy caused by friction of the Earth surface 
and the boundary layer has a considerable value, which is enough to cover the known deficiency 
of the total energy balance of the atmosphere, it becomes clear that in the framework of the 
axion mechanism of the solar dynamo - geodynamo connection it is the solar magnetic field 
that is a host power pacemaker of the Earth global climate [118]. 

7.4 Axion-like particle and extragalactic background light 

A complete physical "portrait" of the axion, obtained from the axion mechanism of Sun lu- 
minosity and the solar dynamo - geodynamo connection, actually predefined the necessary 
conditions of its detection in astrophysics. For example, for the search of the axion with the 
mass rua = 17 eV in radiative decays of the o — )• 77 type, when the decay photons are emitted 
at or near a peak wavelength 

2/ic 

\a = ^ = I459A, (131) 

we used the known observational data obtained during the study of the diffuse extragalactic 
background which, according to [121], are more suitable than the flux from any particular region 
of the sky. 

Surprisingly, our theoretical fit of the spectral intensity of the background radiation (Fig. 20) 
produced by axion decays describes the experimental data in the near ultraviolet and optical 
bands including the data from several ground- and satellite-based telescope observations [175- 
185] with very good accuracy. 

In this sense we may say that our result, obtained within the self-consistent model of Sun 
luminosity and solar dynamo - geodynamo connection, has a strong and evident "experimental" 
support represented by the astrophysical solar-like axions, provided that this is not just a 
fortuitous coincidence of the data (Fig. 20), and the axions with the mass rUa = 17 eV really 
exist. 

7.5 Plausible dark matter candidate: hadronic axion or 
axion-like arhion? 

Along with the identification of the solar- like axion with the 17 eV mass, we also found the 
possible relic axion-like archion with the 2 eV mass in near-infrared band of the extragalactic 
background (Fig. 22), which behaves similarly to the hardronic axion with highly suppressed 
interaction with leptons under certain conditions. 

Discovering both the hadronic axion and the axion-like archion (with similar properties, but 
different physical nature) made us think about what we actually see and identify. Aren't these 
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two particles just axion-like archions in reality? All that we know about the World, tells us that 
Mother Nature, although is notable for its Darwinian diversity in the elementary particles zoo, 
would not "multiply entities beyond necessity". Keeping in mind this principle of Ockham's 
razor and our intuitive notion about the laws of astroparticle physics, let us give some arguments 
in favour of the statement that we are really dealing with the axion-like archions and not with 
hadronic axions. 

The main reason of our doubts is the following. It is known [215] that the model of an 
archion automatically avoids the serious problems related to the necessity of a stable relic 
supermassive quark existence, predicted by the hadronic axion model. As a result, the so-called 
"wild isotopes" should exist as a consequence of a corresponding nucleosynthesis process, and 
they are not observed today, as is known. On the other hand, including the hadronic axion model 
into the Grand Unified model should lead to an inevitable existence of a supermassive lepton 
associated with the axion. This is due to the fact that within such unification, a composite of a 
supermassive quark with an ordinary light quark would cause a supermassive quark instability 
as well as the axion interaction with leptons, which is incompatible with the hadronic axion 
properties [216]. In other words, the axion just would not be hadronic in this case. 

Everything is completely different within the archion model. In contrast to the hadronic 
axion model, the archion model (i.e. the Berezhiani-Saharov-Khlopov (BSK) model of horizontal 
unification [215]) naturally gives rise to a supermassive quark instability and a strong suppression 
of the axion-like archion interaction with leptons [209-216]. At the same time, because of the 
suppression of the archion interaction with the lightest generation of leptons (tt,d,e,z^e), the 
existing limitations on the corresponding scale (the energy scale fpq associated with the break- 
down of the U(l) PQ-symmetry) of the invisible axion are reduced to 

fa = fpQ-W^ GeV, (132) 

which makes the BSK-model of an archion rather close to the hadronic axion model [32-35, 107, 
108]. 

It is interesting that for such values of the energy /„ in the BSK-model of horizontal uni- 
fication the stable massive neutrinos dominance becomes possible, or more strictly speaking, 
the stable sterile neutrinos with the mass about several keV and small mixing angles with the 
active neutrinos. If we also take into account that one or more of these gauge-singlet fermions 
can have Majorana masses below the electroweak scale, in which case they appear as sterile 
neutrinos in the low-energy theory [224], it becomes clear that these particles have fundamental 
importance extending the Standard Model by gauge singlet fermions, which can accommodate 
the neutrino masses. Furthermore, if one of the sterile neutrinos has mass in the 1-20 keV range 
and has small mixing angles with the active neutrinos, such particle is a plausible candidate 
for dark matter [225]. The same particle could be produced in the supernova explosion, and 
its emission from a cooling neutron star could explain the pulsar kicks, facilitate core collapse 
supernova explosions, and affect the formation of the first stars and black holes (e.g., [224] 
and Refs. therein). Therefore, there is a strong motivation to search for signatures of sterile 
neutrinos in this mass range, especially as the first results were obtained recently indicating the 
detection of the sterile neutrino mass of 5.0ib0.2 keV and a mixing angle in a narrow range for 
which neutrino oscillations can produce all of the dark matter [226] . 

We tried to find a trace of the radiative decay of the sterile neutrino in the form i^^ — t- i/ -|- 7 
as well. As we noted before, one should not be embarrassed by the fact that the corresponding 
life-time of the sterile neutrino is many orders of magnitude larger than the age of the universe 
(see (95) and (97)), because sterile neutrinos are produced in the early universe by neutrino 
oscillations [225] and, possibly, by other mechanisms as well, and therefore every dark matter 
halo should contain some fraction of these particles. In other words, given a large number of 
particles in these astrophysical systems, even a small decay width can make them observable 
via the photons produced in the radiative decay. 

We limited our search for possible 7-spectra of the sterile neutrino decay to the analysis of 
the experimental results obtained by different authors during their study of the diffuse X-ray 
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Figure 25: (a) The cosmic XRB spectrum and the predicted contribution from the active galactic nuclei (AGNs), 
that gives origin to the X-ray background. Gray points: HEAO-l A2 HED data (Gruber et al. [227]). Dark green 
points: HEAO-l A4 LED (Gruber et al. [227]). Cyan points: Rossi-XTE (Revnivtsev et al. [228]). Red bowtie: 
HEAO-l A4 MED (Kinzer et al. [229]). Blue bowtie: ROSAT PSPC data (Georgantopoulos et al. [230]). Light green 
bowtie: BeppoSAX (Vecchi et al. [231]). Purple and yellow bowties: Newton- XMM (Lumb et al. [232]; De Luca & 
Molendi [233] ) . Solid line: synthesis model spectrum by Haardt-Madau [234] , produced by a mixture of absorbed (the 
short-dashed black line) and unabsorbed (the long-dashed black line) AGNs. Adopted from [234]. 
(b) The cosmic XRB spectrum as a function of the wavelength Aq. Solid purple line: synthesis model spectrum by 
Haardt-Madau [234], produced by a mixture of absorbed and unabsorbed AGNs. Black points with bars (crosses): 
HEAO A-2 results (Wu et al. [235] ) . Solid green line: our fit of cosmic XRB spectrum. The orange region is a sum 
of our fit of spectrum (the green line) and the spectrum of 7-quanta (the blue shaded region) born in the supposed 
radiative decay of the sterile neutrinos. 

background (Fig. 25a). The cosmic XRB spectrum in the 0.1-1 keV range was particularly 
interesting in this sense, because the X-ray spectrum here is not "distorted" by the absorption 
processes, since unabsorbed AGNs dominate in this energy region (see predictions by the Haardt- 
Madau model on Fig. 25a). Allowing for the fact that in this energy region the problem of 
Lyman-a forest may be neglected, it becomes possible to use Eq. (79) for theoretical calculation 
of the intensity of sterile neutrino contributions to the diffuse X-ray background as a function 
of the wavelength Aq. 

It is necessary to mention some specifics of cosmic XRB spectrum calculation for this re- 
gion^^. First, the XRB spectrum is highly "damaged" (see Fig. 25a) by multiplicity of observa- 
tions with different measuring bases, and possibly, different methodologies. For this reason we 
decided to use the observational data obtained in the single HEAO experiment. It covers the 
whole spectrum in 1 - 400 keV range, and also provides the measurements of the cosmic XRB 
spectrum in 0.1 - 1 keV region [235]. Second, when calculating the 7-spectrum formed by the 
radiative decay of sterile neutrinos, we used Eqs. (79)-(81) adjusted for expressions (95)-(97), 
which reflect the radiative features of sterile neutrinos, and assumed their velocity distribution 
(see (81)) to be Gaussian with a variance of Vc ~ 30 km/ s [236, 237] and their density to be 
equal [238, 239] 



^shl 



0.3 



sin^ 29 



ms 



(133) 



10-10 ) VlOO keV. 
We also introduced a factor of 3 for the cosmic XRB spectrum calculation (79), which 



i^The observations performed using the X-ray spectroscopy of locations in the Universe are not used here, because 
we calculate the spectrum of 7-quanta born in radiative decays 1/5 — ?> z/ -h 7 for different z (integral method (86) [172]) 
in contrast to the X-ray spectroscopy of emission lines from the cooling stars (so-called differential method (e.g. [162]) 
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manually takes into account that the bulk of the EBL contributions in the decaying-neutrino 
scenario conies from neutrinos which are distributed on larger scales. We will refer to these 
collectively as free-streaming neutrinos, though some of them may actually be associated with 
more massive systems such as clusters of galaxies, and a possible impact of the astrophysical 
data non-gaussianity. In support of this assumption let us quote a remark made in paper [240], 
where a factor of 10 was used: "...the fact that a significant part of the dark matter at redshifts 
z ^ 10 is concentrated in galaxies and clusters of galaxies just means that the strongest signal 
from the dark matter decay should come from the sum of the signals from the compact sources 
at z ^ 10. Taking into account that the DM decay signal from z ^ 10 is some two orders of 
magnitude stronger than that from z ^ 10, while the subtraction of resolved sources reduces 
the residual X-ray background maximum by a factor of 10, we argue that it would be wrong to 
subtract the contribution from the resolved sources from the XRB observations when looking 
for the DM decay signal" . 

Our theoretical fit (the yellow curve in Fig. 25b) of spectral intensity of the background 
radiation (79) produced by sterile neutrino decays describes the experimental data from HEAO 
A-2 in near X-ray band [227] adequately. However, this result requires a serious and thorough 
verification, since the data by Wu et al. [235] are, in fact, the only cosmic XRB spectrum 
measurements of the Large Magellanic Clouds in the 0.16 - 3.5 keV region. From this point of 
view, this result is more likely to serve as a demonstration of the possibilities and peculiarities 
of the sterile neutrino search in the diffuse extragalactic^*^ X-ray background which sometimes 
may be more suitable than the flux (in the form of a narrow band) from any particular region 
of the sky. 

Turning back to the archion properties, let us point out that such sterile neutrino {ms = 
2.75 keV, sin^ 26 = 5.1-10"®) representing the dark matter with density Og ~ 0.20 (see (133) and 
Fig. 26) fits into the archion model very well. This is due to the fact that it is the sterile neutrino 
appearance that makes it possible to solve the problem of missing dark matter (see (92)) which 
is associated with the archions model (in a case of ms > nia) and emerges in LTR cosmology: 

^DM = ^s + K"" ~ 0.25, (134) 

where OJJ*'^ ~ 0.04 is a non-standard density of the axion dark matter (92). 

Briefly summarizing all said above about the problem of hadronic axion and axion-like 
archion, in our opinion, the archion model is preferable as compared to the hadronic axion 
model, but as the saying is, time will tell. 

And finally let us emphasize two the most painful points of the present paper, at least from 
its authors' point of view. 

Regardless of the model type, the simultaneous existence of two axion-like particles (the 
hadronic axion with the mass 17 eV and the axion-like arhion with the mass 2 eV) rises a natural 
question whether there may be two different energy scales fa associated with the breakdown of 
the U(l) PQ-symmetry, one of them coming from the hadronic axion (fa = 3.5 • 10^ GeV~^), 
and another relic one - from the axion-like archion (fa = 3.0 • 10^ GeV~^). Or otherwise 
stated, can the energy have two scales, or in general, be hierarchical? If it can, then what kind 
of consequences of such fundamental symmetry violations hierarchy may be observed today? 
However, if one just winks at the possible existence of the axion-like particle with the 2 eV 
mass, the problem just vanishes without any consequences for the results of the present paper. 

The second painful point is related to the key problem of the axion mechanism of Sun 
luminosity and is stated rather simply: "Is the process of axion conversion into 7-quanta by the 
Primakoff effect really possible in the Solar tachocline magnetic field?" This question is directly 
connected to the problem of the hollow magnetic flux tubes existence in the convective zone 
of the Sun, which are supposed to connect the tachocline with the photosphere. So, both the 
theory and experiment have to answer the question of whether there are the waveguides in the 
form of the hollow magnetic flux tubes in the convective zone of the Sun, which are perfectly 

Wu's analysis [235] of the background data reveals a limit on the mean absolute intensity for the extragalactic 

-8 „^„„ „^,-2„-l„^-l 
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Figure 26: Full parameter space constraints for the sterile neutrino production model based on the diffuse extragalactic 
X-ray background, assuming sterile neutrinos constitute the dark matter. To facilitate comparisons, we adopt many of 
the conventions used by Abazajian [241]. Favored regions are in red/magenta colors, disfavored and excluded regions 
are in blue/turquoise colors. The favored parameters consistent with pulsar kick generation are in horizontal hatching 
(Kusenko et al. [242, 243]). Constraints from X-ray observations include the diffuse X-ray background (turquoise) 
(Boyarsky et al. [240]). Also shown is the best current constraint from Chandra, from observations of contributions of 
dark matter X-ray decay in the cosmic X-ray background through the CDFN and CDFS (brown contour, "Unresolved 
CXB Milky Way"). Also shown is an estimate of the sensitivity of a 1 Ms observation of M31 with 1X0 (yellow). The 
region at ms < 1.7 keV is disfavored by conservative application of constraints from the Lyman-a forest. 
Inset: Supernova bound on sterile neutrino masses ms and mixing angles sin (20), where the purple region is excluded 
by the energy-loss argument while the green one by the energy-transfer argument [238]. The excluded region will be 
extended to the dashed (red) line if the build-up of degeneracy parameter is ignored, i.e., r/(t) — 0. The dot-dashed 
(green) line represents the sterile neutrinos as dark matter with the correct relic abundance ^s^q — 0.1. The red star 
{ms = 2.75 keV, sin^ 29 = 5.1 • 10~®) marks our result of sterile neutrino parameters identification. 

transparent for 7-quanta, or our model of the axion mechanism of Sun luminosity and solar 
dynamo - geodynamo connection is built around simply guessed rules of calculation which do 
not reflect any real nature of things. 
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A Appendix I. Effect of 7-quanta channeling in 
periodic structures 

It is known [1] that the real part of dielectric susceptibility xi^) for the photons with energies 
exceeding the X-electrons binding energy has the form 

Rexiuj) = -WeV^', (A.l) 

where coe = (47rA^e^/me)^" is the electron plasma frequency, N is the electrons density, rue is 
the electron mass. Since cOg ~ 10 eV for the majority of materials, the susceptibility is small and 
negative in the X-ray band. It means that X-ray photons may experience the total reflection 
on the border of two materials from the one with the larger value of \Re x{^)\i ^^^ the angle 
of photons arrival to the border should not exceed 9c = \Re x|^ > where A{Re x) is a dielectric 
susceptibility step. It is appropriate to mention that the phenomenon of a small-angle X-ray 
reflection has been used in X-ray optic elements for a long time [2], and also for transporting 
(e.g. [3]) and turning (e.g. [4]) the X-ray bundles by cylindrical tubes. 

On the other hand, it is also known that according to the optical-mechanical analogy (Hamil- 
ton and Fermat principles identity) , in certain cases the radiation propagation may be described 
in terms of ray trajectories obeying the principle of least action. The refractive index, defining 
a ray trajectory, plays a role of an external potential in this case. Based on this analogy, one 
may conjecture that there exists an effect of high-energy chargeless particles channeling that 
leads to the substantial anisotropy after passing the periodic structures. This assumption is 
validated using the example Vinecky and Finegold model problem [5, 6] related to calculation 
of the ray trajectories and absorption coefficients for the hard photons in geometrical optics 
approximation. 

A.l Statement of a problem 

Let us follow [5, 6] and consider 7-radiation with frequency uj = 2ttc/X transmission through a 
medium under the following condition: 

Ae < A < ao, (A.2) 

where Ag = h/mc is the electron Compton wavelength, ao is a typical interatomic distance in 
the medium. 

The left-hand side of the inequality (A.2) lets one consider the quanta as propagating in the 
continuous medium with a refraction index n. As long as the right-hand side of (A.2) holds, the 
incident wave frequency is substantially higher than the lattice vibration eigenfrequencies (at 
least for the higher electron shells of the atoms that form the lattice) . Therefore the scattering 
electrons may be considered free, and a crystal in the large may be treated as a frozen spatially- 
inhomogeneous electron plasma (the impact of nuclei on the scattering is negligible) . 

The electron plasma dielectric constant^^ may be written down in the form 

e = ei+ie2 = l- .'^'^''^ , (A.3) 



^^For a justification of the term "dielectric constant" applicability for the X-ray and 7-bands in the absence of the 
Lorentz field averaging over the physically small volume elements see [7] . 
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Figure A.l: Ray trajectory dependence on the angle a in a long-period structure: the channeling happens when 
a < a„j; when a > «„, the ray crosses the surfaces. 

where w-y is the damping parameter which defines the wave absorption in the system. 
For the refraction index 

n = "v/i = 77-1 + in2 
from (A. 3) taking into account that oj ^ Ue ^ oj-y in the 7-band, we derive 



ni ~ 1 



1 /w.\2 



n2 



1 Uj'^iO'y 



(A.4) 



2 Vw/ ' "^ 2 a;3 ■ 

In order to perform the quantitative calculations, let us consider a simplified model of the 
one-dimensional lattice (a stack of alternating layers with different, but constant, densities and 
uniform N(r) distribution along the layer (Fig. A.l)). 

Let us choose the distribution in the y direction (perpendicular to the layers) in the form: 



N(y) = No(l + ^^sm^TT- 
V a 

Corresponding expressions for ni, 71-2 are 



(A.5) 



m = 1 - '^e (1 + /?' sin' l^y, ns = ^CC' (l + /3' sin^ h^y 



where 



(A.6) 






Wn A-KNoe^ 



OJ^ 



<1, C= — ^i. ka = 



27r 



(A.7) 



So the problem reduces to determining the intensity of the 7-quanta beam hitting a sample 
of the thickness x with the angle a to the Ox axis (Fig. A.l). 

We shall use the equation for the intensity of a ray that passed a path / in the absorbing 
medium: 



Ji = Joexp < —2k / n2{l)dl > , k 



27rni 27r 



(A.8) 



The integral in (A.8) is taken along the ray trajectory, while the trajectory is determined 
by means of the Fermat principle [8] 

5 / nidi = 0, 



variation of which yields a system of three nonlinear second-order differential equations describ- 
ing the ray equation in a parametric form x = x(l), y = y{l), z = z{l): 



ni- 



d'^Xi 
dp 



+ 



drii dxj\ dxi dni 



dxj dl 



dl 
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dxi 



i,j = 1,2,3 



(A.9) 



(xi, X2, X3 correspond to x, y, z). It follows from the identical relation dxjdxj = dt^ that only- 
two equations of (A. 9) are independent. 



A. 2 Solution of the model problem 

For the model (A. 5) the ray trajectory lies in one plane xy, and the system (A. 9) reduces to 
the equation 



d'^y dni /dy 
""^ dp dy \dl 



drii 
dy ' 



(A.IO) 



Let us point out some features of the channeling process that directly follow from Eq. (A.IO). 
Assuming dy/dl = everywhere along the ray, we derive the particular solution (A.IO) which 
describes the propagation in the planes parallel to the planes of constant density. According 
to (A.IO), such motion is only possible when dni/dy = 0, i.e. within the planes Aj and Bj^ 
corresponding to the minimum and maximum values of ni(y). Otherwise, the straight motion 
of a ray along the planes ni{y) = const is impossible. Particularly, if the ray is initially parallel 
to the planes, but lies neither in Aj, nor in Bj, it is bent towards the nearest Bj, since ni 
grows in this direction. Therefore, the motion within the plane Bj is stable, while within Aj it 
is not. In other words, the rays entering the layers along the Ox axis are pushed out into the 
"channels" of the low electron density (regions with minimum N(y)). 

The first integral of Eq. (A.IO) is as follows: 



dy 
dl 



1 



ni{yo] 

ni{y) 



cos a 



1/2 



c?y 
dx 



ni{y) 



ni{yo) cos a 



1 



1/2 



(A.ll) 



Let us use (A.ll) for determining the area of the ray motion depending on initial values of 
a and yo- By taking dy/dx = in the ray's turning point y = ym in (A.ll), we write down 



nliVm) = ni(yo)cos^a. 



By substituting (A. 6) we derive 



1 



„ sin^ a + cos^ a ■ sin^ ( -kayo ] ■ 



sin -kay„i , — ,^^n - . --- -- , 

Eq. (A. 13), obviously, has a real root ym only in the case 



1-e 



1 



sin a + cos a • sin ( -kayo ] ^ 1, 
whence taking into account (A. 7) we obtain 

a ^ am(yo), am(yo) = arctan 
The dy/dx function zeros have the form 

ym = ±A;^^ arccos (cos kayo - 2(/30 



/3^cos ( -kayo 



<1. 



-2 • 2 

sm a 



)■ 



(A.12) 



(A.13) 



(A.14) 



(A.15) 



(A.16) 



For each yo a ray with the angle of arrival a < amivo) "oscillates" between the end points 
[yrm—ym] when passing through a sample, i.e. it "channels" between two neighbouring planes 
Aj and Aj^i. Hence, Omivo) is the maximum angle of channeling for the given yo- According 
to (A.15), the angle Um = corresponds to yo — Vm = (l/2)a, which means that the values of 
yo that make the channeling possible lie in the region 

11 

-a, -a 
2 '2 
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thus embracing the whole sample facet exposed to the beam. 

If a > amino), the dy/dx function has no zeros ~ a ray that entered under the angle 
a > amino) does not channel - instead it crosses the layers one by one and is heavily absorbed 
by them. Fig. A.l shows the trajectories of both the "captured" ray that propagates between 
two planes throughout the sample and the non-channeling one. 

From this analysis the following picture of the channeling emerges. When a beam hits the 
sample in the q = direction, virtually all the rays pass through the sample without any 
substantial absorption - including the ones that propagate within the strongly absorbing layers 
as they are pushed out back into the channels. Under small a ^ the thin bands of 5yo values 
appear on both sides of the low density planes, where the incident rays quit channeling by 
switching to the crossing trajectories. Since each crossing is accompanied by a loss of some 
amount of energy, we obtain two sets of rays - experiencing the weak and strong absorption 
respectively. With a growth (and, consequently, the 6yo bands growth) first set of rays shrinks, 
while the second set grows. When a = amivo = 0), the channeling region collapses and all the 
rays are strongly absorbed. Therefore the value 

um = amivo = 0) = arctan/3^ w /3^ — /3 — (A. 17) 

characterizes the maximum possible channeling angle, thus determining the divergence of a 
beam for a given sample. As seen from (A. 17), the value of aM decreases with the beam 
rigidity growth! 

Let us estimate the value of aM for the plasma medium in the Solar tachocline zone. Let 
/? ~ 1 for simplicity. Taking into account that the plasma frequency in the tachocline is equal 
to u)e ~ 4.6 • 10^^ s~^ [see (8) in the main part], and the frequency of 7-quanta with the 
energy {E) = 4.2 keV (or A ~ 4.6 • 10"^ cm) is w = wttc/\ = 2TT{E)/h ~ 4 • lO^^ s'^ , from 
(A. 17) we obtain am ~ 10"^. For the sake of channeling effect illustration. Fig. A.l shows the 
model results of the photon (A ~ 10~^ cm) beam propagation through the layered media with 
aM ~ 10~^ (see Section A. 4). 

Let us now integrate the Eqs. (A. 11). For the case of (A. 4) these equations are reducible to 
the following form (up to the first vanishing terms of the order ^^ under the radical sign) : 

dy . L 2 ■ 2 f^r \ dy 2 ■ 2 A 

di y 

where 



psinaWl-g^sin^/ fcayj, — = ptanaWl - g^ sin^ ( -fc^y ) , (A.18) 



/-, ,r>^^n o .9/^1, \ /3£cota /. o TTWo tan^a\ /. ,„n 

p=s\l + iPi)^coi^asin^{-kayoy q=— =(sm2^ + -— ^j . (A.19) 

By integrating (A.18) we obtain the following trajectory equation: 

^^^ I 2 /F(arcsingsin(^fcay),g-i) at g ^ 1, 
' kaPi\qFi\kay,q) at q < 1 

Here Fiip, q) is an elliptic integral of the first kind, and the inverse function amF (the Jaco- 
bian elliptic function) period is equal to 4i^(g), where Kiq) = FiTT/2,q) [9]. Correspondingly, 
the function y(x) inverse to (A. 20a) is periodical with the period 

., = ^^^-(^-') (A.21) 

and is associated with the channeling trajectories. For the crossing trajectories (A. 29b) which 
are translation-invariant under the simultaneous transformations 

x-^x + jxr, y-^y + ja, j = 0,1,2,..., (A.22) 
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we have 






KPi 



K{q). 



(A.23) 



The regions of channehng and crossing trajectories are dehmited by the value of the angle 
a = ctmiuo) corresponding to the value q = 1. When a = am{yo), the rays asymptotically 
approach the "repelling" planes Aj and x,- — )• oo. 

A number of trajectory oscillations along the distance x for the channeling and crossing rays 
is, respectively 



Nr 



X 

Xr 



1 U-(.-;) atq>l, 

8 "^^{q-^K-^q) at q<l. 



(A.24) 



A. 3 Determination of the absorption coefficient angular 
dependence 

Let us move on to the absorption coefficient calculation in the exponent (A. 8). By means of 
(A. 18), integration over trajectory reduces to integration over a variable u = kay/2 within the 
region 

{0 ^ n ^ Um at g > 1, 
^ n ^ 7r/2 at q<l. 

Let us remind that u = kaVm/'^ = aicsui{q~^) by virtue of (A. 15) and (A. 19) corresponds 
to the turning points ym of the channeling trajectories. Further, substituting (A. 8) into the 
expression for n2 (A. 6), performing integration and multiplying the obtained integrals, which 
correspond to the above mentioned regions, by A^^^ (A.24), we obtain 



O" 




E(q) ' 
Kiq) 



at q > 1, 
at q < 1, 



(A.25) 



where E{q) is the complete elliptic integral of the second kind [10]. 

According to (A.25), the intensity J{x) (see (A. 8)) of a ray that passed through a sample of 
the thickness x may be written down as 



J(x) = Jo exp(— o"x) = Jo exp 



Xo 



cos a 



x) ■Qia,yo,x), 



(A.26) 



where 



Q{a,yo,x) 



exp 



exp 




q^ cos a 



at q > I, 
at q < 1. 



(A.27) 



Here the multiplier Joexp(— xoa^^/cosa) corresponds to 7-rays propagation through a homo- 
geneous medium with the electron density A^e and the absorption coefficient xo = ^C'C^- An 
additional multiplier Q{a,yQ,x) characterizes the influence of the medium layering. Fig. A. 2 
shows the Q{a,yo,x) curves for 



= 0, 1/8, 1/4, 3/8, 1/2 (A.28) 

2, obtained by a numerical calculation using the elliptic 



m 

a 

with A = 10^^ cm, /3 = 1 and xo^ = 
integral value tables [10]. 

As seen in the figure, the "total transmission" (Q = 1) is observed when a = 0, yo = 
(the propagation in the Bj plane), and the maximum absorption depending on f3 and {Q — )• 
exp(— xo/5a^cosa)) is observed when a = 0, yo = (l/2)« (the propagation in the Aj plane). 
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Figure A. 2: Transmitted beam relative intensity dependence on the arrival angle a for different values of j/o (blue 
curves; the symmetric parts of the curves for a < are not shown). Red curve represents Q{a,yQ,x) averaged over 
all values of j/o • 

The dips on the decaying parts of the transmission curves when yo / (l/2)a correspond to 
the critical values of the angle a = am{yo) (A. 15), for which the 7-quanta beam, that entered a 
sample, asymptotically approaches the nearest Aj plane and is absorbed in it. With somewhat 
bigger values of a the rays pass on to the neighboring inter-plane space and moving within it 
leave a sample before they enter the next ^j+i plane. The additional narrow maxima on blue 
curves (Fig. A. 2) correspond to this type of rays. 

Finally, when a > aM and, particularly, when a — )• 7r/2 (propagation across the layers) 
from (A. 18) we obtain g — )• 0, which corresponds to Q — )• ex.p{—xo(3'^x/2cosa), i.e. there is an 
additional absorption due to the averaged additional density of the layers A''! = Nq/3'^. 

Fig. A. 2 shows a curve Q{a, yo, x) averaged over all values of yo- The side dips and maxima 
are smoothed out, and there is only a central maximum. It corresponds to a clearly defined 
primary 7-quanta propagation with their angle of arrival within ajv/ ~ 2 • 10~^. Therefore, 
the dense layers "modulating" a material play a role of an anisotropic filter which passes the 
7-radiation for a narrow range of angles a < aM only. The physical mechanism of the arising 
transparency anisotropy consists in 7-quanta channeling between the layers because of the 
radiation "refraction" in a heterogeneous electron plasma. 

A. 4 Channeling effect onset conditions 

By definition, channeling occurs for the particles the motion of which in the transversal phase 
plane is limited to a region 



AyAp ~ aahk, 



(A.29) 



where a is a channeling angle. According to the uncertainty principle, the size of this region 
cannot be less than h, i.e. 



akn ^ a 



(A.30) 
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By letting a ~ We/w for the estimation, from (A. 30) we derive 

ako ^1, ko = ujq/c. (A. 31) 

For example, in a monocrystal with a = oq ~ 3 • 10~^ cm, cjq ~ 5 • 10^^ s~^ we have 
aoko ~ 5 • 10~^ <C 1, i.e. the condition (A. 31) does not hold and, consequently, the channeling 
is impossible. This prohibition is true for any quanta, since in the region a; ^ cjq, where the 
approximation (A. 4) and (A. 5) is applicable, the frequency uj falls out from the condition (A. 31). 

However, for a layered system with a » oq the condition (A. 31) may prove to be true, 
and the channeling may be possible. A more strict condition ako ^ 1 makes the classical 
description of this effect possible. This is the case for a long-period structure with a/ao S> 10^ 
considered above. It is appropriate to emphasize that we consider here the structures based 
on the amorphous matrices which are not monocrystals. Otherwise the long-period structure 
would have been overlapped by the short-period oscillations with the amplitude equal or greater 
than PNq, which would have lead to a noticeable tunnel effect. The question about channeling 
possibilities in this case requires a special research. 

Thus, the papers [5, 6] show that the phenomenon of X-ray and 7-radiation channeling exists 
in layered structures under conditions when it is possible to use geometric optics. The essence 
of this phenomenon lies in the fact that the rays are reflected form the layers with the higher 
electron density if they propagate under the small enough angle (a < om) to the layers plane. 
The initially uniform intensity distribution in the transversal plane becomes non-uniform, since 
the rays concentrate in the "channels" - the layers with the lower electron density. It leads to 
the substantial absorption decrease and deeper radiation penetration into the sample along the 
layers than in the case of an arbitrary arrival angle. 

A. 5 On the account of an absorption impact on X-ray 

intensity when channeUng through the solar layered 
structures 

As it was shown above, the process of 7-rays channeling through a homogeneous medium with 
the electron density N^ and the absorption coefficient xo = ^C'^^ may be described by the 
multiplier Joexp(— xoa^/cosa) in (A. 26). 

Calculation of this multiplier for an arbitrary point in the solar convective zone, obviously, 
requires the estimation of the average Rosseland free path or so-called Rosseland opacity (the 
photon absorption coefficient averaged according to Rosseland) in these points. On the other 
hand, it is necessary to know the radial profiles of the temperature and density in the solar 
convective zone in order to calculate the Rosseland free path or Rosseland opacity. 

The transmission of photons with the intensity Jo normally incident on a uniform plasma (A. 26) 
is given by 

T{u) = J{u)/Jq{v) = exp(-xoa;) = exp [-k{u)px\ , (A.32) 

where hv is the photon energy and J{v) is the attenuated photon intensity emerging from the 
plasma, k{v) is the opacity per unit mass typically measured in units of err? jg, p is the density, 
and X is the optical path length. For plasmas such as the Sun that are much larger than the 
photon mean free path, radiation transport is usually described by the diffusion approximation 
[11-13] using the Rosseland mean opacity kn, 

^T3 I r ^1T3 

(A.33) 

where B is the Planck function, T is the plasma temperature, and the weighting function dB/dT 
peaks at roughly 3.8 kT. Near the convective zone (CZ) boundary T ~ 190 eV and dB/dT 
peaks at hu ~ 750 eV (Fig. A. 3). The frequency dependent opacity near the CZ boundary 



76 




1000 




I i_ I :. I I I I I I I I I I I I 

B = Planck fiincion . 
^T = temperature 



opacity "window" 
I I I I I I I I I I I I 




JLJ 



1.0 



H0.8 

0.6 dB 
dT 
H0.4 



500 1000 1500 2000 2500 3000 
hv (eV) 



0.2 
0.0 



Figure A. 3: Frequency dependent opacity [14-16] for the 17 element solar composition [17] near the base of the 
solar convection zone compared to dB/dT. The electron temperature and density were 193 eV and 1 • lO^'^ cm~^, 
respectively. Adopted from [18]. 

calculated using the opacity project model [14-16] is displayed in Fig. A. 3. Comparison with 
the weighting function for the Rosseland mean shows that the most important photon energies 
are approximately 300 < hv < 1300 eV. 

However, it is necessary to point out some essential peculiarities of the absorption impact on 
the X-ray intensity when it channels inside the solar layered structures based on, e.g. magnetic 
flux tubes superlattices (see Appendix B). 

First , as it was noted in the body text, the total energy balance of the Sun is not violated in 
the framework of the axion mechanism of Sun luminosity, but the radiation transport changes 
substantially relative to the standard model of the Sun - the part of radiation transport not 
related to axions is very small (~ O.OISA^^^). It means that the thermodynamic parameters 
(temperature, pressure, plasma density, electron density etc.) are considerably smaller in the 
axion model of the Sun as compared to the standard model. On the other hand, the Rosseland 
free path or Rosseland opacity calculation inside the thin magnetic flux tubes, which play the 
role of 7-quanta waveguides formed in the tachocline, must take into account the new 
values for the mentioned parameters in the framework of the axion mechanism of Sun 
luminosity. 

Second , we suppose that the decrease of the pressure in the convective zone, for example, by 
an order of magnitude may lead to such decrease of the pressure inside the magnetic flux tubes 
that it virtually does not influence the radiation transport in these tubes. In other words, the 
Rosseland free paths of 7-quanta in the thin magnetic flux tubes is so big that the corresponding 
Rosseland opacity tends to zero, and so do the absorption coefficients in (A. 32). The low 
refractivity, or equally the high transparency of the thin magnetic flux tubes is achieved due to 
the high magnetic pressure (see (16) and Fig. 9) able to compensate the outer pressure of the 
convective zone completely (see Appendix B for details). 
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B Appendix II. On a possibility of the layered 

structures formation in the solar convective zone 
on the basis of the magnetic flux tubes 
superlattices 

It is natural to ask a question, whether there is a physical possibility for the above mentioned 
long-period structures formation in the magnetohydrodynamical plasma media typical for the 
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solar dynamo evolution. Curiously enough, such mechanisms do exist. Let us examine some of 
them briefly. 

B.l Zonal jet streams 

The long-period structures in magnetohydrodynamical plasma media may show as the so-called 
zonal jet streams, spontaneously generated in turbulent systems. In fact, zonal jets are very 
common in nature. Well-known examples are those in the atmospheres of giant planets and the 
alternating jet streams found in the Earth's world ocean [1]. Zonal flow formation in nuclear 
fusion devices are also well studied [2]. 

As we have already pointed out above, a common feature of these zonal flows is that they 
are spontaneously generated in turbulent systems. Because the Earth's outer core is believed 
to be in a turbulent state, it is possible that there is a zonal flow in the liquid iron of the outer 
core. It is interesting that a previously unknown convective regime of the outer core that has 
a dual structure comprising inner, sheet-like radial plumes and an outer, westward cylindrical 
zonal flow^^, was recently found [4] (Fig. 6 in the body text). 

Fig. 6 in the main body of the present paper shows snapshots of the same data of zonal flow 
formation in the Earth's core. The sheet convection structure is visualized as isosurfaces of the 
axial vorticity, which are almost straight in the z direction. The blue curves surrounding the 
iCz isosurfaces are streamlines of the velocity in the outer part of the dual-convection structure. 
The ring- like shape of each streamline indicates that the azimuthal component is dominant. 

However, according to (A. 26) in A, for the ideal (without absorption) photon channeling 
the long-period structures, in which one of the interlaced media has almost zero density, are 
necessary. Surprisingly, such long-period (in terms of density) structures may appear in the 
plasma media in general, and in the solar convective zone in particular. We mean here the 
so-called magnetic flux tubes, the properties of which are discussed below. 

B.2 Some properties of the magnetic flux tubes in Sun 
convective zone 

Fig. B.l shows the results of the three-dimensional solar hexagonal magnetoconvection simula- 
tion [5] . In addition to a mere fact of the long-period layers formation based on the magnetic 
flux tubes superlattices, it is necessary to make sure that the matter density inside the tubes is 
much smaller than the density of the outer plasma. 

At the same time, as the analysis in Appendix A shows, for the purposes of photon channeling 
it is not necessary to guarantee the layers periodicity. In other words, the channeling requires 
a large number of necessarily interlaced, although not necessarily periodic, layers of different 
density. The total crosscut area of such interlaced layers must also be comparable to the photon 
beam cross-section. 

In this connection below we shall examine some properties of the magnetic flux tubes in the 
Sun convective zone, which may serve as a basis for estimating the temperature, pressure and 
matter density inside the tubes depending on the similar plasma parameters outside the tube. 

^^ Computer simulations have been playing an important role in the development of our understanding of the 
geodynamo, but the direct numerical simulation of the geodynamo with a realistic parameter regime is still beyond 
the power of today's supercomputers. Difficulties in simulating the geodynamo arise from the extreme conditions of 
the core, which are characterized by very large and very small values of the non-dimensional parameters of the system. 
Along them, the Ekman number, E, has been adopted as a barometer of the distance of simulations from real core 
conditions, in which E is of the order of 10^^^. Following the initial computer simulations of the geodynamo, the 
Ekman number achieved has been steadily decreasing, with recent geodynamo simulations performed with E of the 
order of 10^^ [3]. In work by Miyagoshi et al. [3, 4] they present a geodynamo simulation with the Ekman number of 
the order of 10^^ - the highest resolution yet achieved, making use of 4096 processors of the Earth Simulator. And 
what is ahead when the magnitude of E becomes closer to its real value?! 
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Figure B.l: Vertical field component for Rm = 400 when the imposed field is vertical, displayed by perspective plots. 
Four times are shown, for the three levels z = 0, 1/2, and 1. In this and subsequent plots the time unit is (cell 
height)/(inax. vertical velocity), i.e. ^ 1.4 turnover time. Adopted from [5]. 

We shall also derive the zero refractivity (or absolute transparency) condition for the effective 
photon channeling along the magnetic tubes. 

B.3 The self-confinement of force-free magnetic fields and 
energy conservation law. 

Magnetic field B alternating along the vertical axis z induces the vortex electric field in the 
magnetic flux tube containing a dense plasma. The charged particles rotation in plasma with 
the angular velocity co leads to a centrifugal force per unit volume 



~ p\uj\ r. 



(B.l) 



If we consider this problem in the rotating noninertial reference frame, such noninertial- 
ity, according to the equivalence principle, is equivalent to "introducing" a radial non- uniform 
gravitational field with the "free fall acceleration" 



9{r) 



\CJ\ r. 



(B.2) 
liquid" of 



In such case the pressure difference inside {pmt) and outside {pext) the rotating 
the tube may be treated by analogy with the regular hydrostatic pressure (Fig. B.l). 

Let us pick a radial "liquid column" inside the tube as it is shown in Fig. B.l. Since the 
"gravity" is non- uniform in this column, it is equivalent to a uniform field with the "free fall 
acceleration" 



{9{r)) 



1 



uJ\'^R. 



(B.3) 



where R is the tube radius which plays a role of the " liquid column height" in our analogy. 

By equating the forces acting on the chosen column similar to the hydrostatic pressure 
(Fig. B.l), we derive: 



Pext =Pint + -/o|wpi?^. 



(B.4) 
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Figure B.2: Representation of the "hydrostatic equihbrium" in the rotating "hquid" of a magnetic flux tube. 

This rises the natural question as to what physics is hidden behind the "centrifugal" pressure. 
In this relation, let us consider the magnetic field energy density 



I-BI 



WB 



2fio ' 



(B.5) 



where ^o is the magnetic permeability of vacuum. 

Suppose that the total magnetic field energy of the "growing" tube grows linearly between 
the tachocline and the photosphere. In this case if the average total energy of the magnetic field 
in the tube transforms into the kinetic energy of tube matter rotation completely, it is easy to 
show that 



T T I R|2 rL-;|2 



(B.6) 



where I = rriE? /2 is the tube's moment of inertia about the rotation axis, m and V are the 
mass and the volume of the tube medium respectively. 
Then from (B.6) it follows that 



p\uJ\R^ 



\B\ 



2 2^0 

Finally, substituting (B.7) into (B.4) we obtain the desired relation 

\B? 

Pext = Pint + 



2fIo ' 



(B.7) 



(B.8) 



which is exactly equal to a well-known expression by Parker [6], describing the so-called self- 
confinement of force- free magnetic fields. 

B.4 Hydrostatic equilibrium and a sharp tube medium cooling 
effect 

Assume that the tube has the length l{t) by the time t. Then its volume is equal to S ■ l{t) and 
the heat capacity is 



c-p{t)-Slit), 



(B.9) 



where S is the tube cross-section, p{t) is the density inside the tube, c is the specific heat 
capacity. 
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If the tube becomes longer by v{t)dt for the time dt, then the magnetic field energy increases 



1 \B\ 



Sv{t)dt. 



(B.IO) 



2 2//0 

where v{t) is the tube propagation speed. 

The matter inside the tube, obviously, has to cool by the temperature dT so that the internal 
energy release maintained the magnetic energy growth. Therefore, the following equality must 
hold: 



dT 

cp{t)l{t) — S 



l\B\ 



2 2/io 



v{t)S. 



(B.ll) 



Taking into account the fact that the tube grows practically linearly [7], i.e. vt = I, the 
Parker relation (B.8) and the tube's equation of state 






p 



M* Pint 



R* T{t) ' 
the equality (B.ll) may be rewritten (by separation of variables) as follows: 



dT 



R^ 

2C;U=f 



Pext 
Pintit) 



dt 

T' 



(B.12) 



where /n* is the tube matter molar mass, i?* is the universal gas constant. 
After integration of (B.12) we obtain 



In 



m 

r(o) 



2c/x=i, 



Pext 

Pintir) 



dr 

T 



(B.13) 



It is easy to see that the multiplier (l/r) in (B.13) assigns the region near r = in the 
integral. The integral converges since 



lim 

r-S>0 



Pext 



1 



0. 



Pintir) 
Expanding pint into Taylor series and taking into account that Pint{T = 0) = Pext, we obtain 

dpintir = 0) 



Pint[T) =Pext + 



dT 



where 
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From (B.14) it follows that 



dT 



Pint{r = 0) 



dT 



Pext 
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1 « 7r, 



Pint{r) 1-JT 

therefore, substituting (B.16) into (B.13), we find its solution in the form 



T{t) = r(0) exp 



2c fl^: 



It 



(B.14) 



(B.15) 



(B.16) 



(B.17) 



It is extremely important to note here that the solution (B.17) points at the remarkable fact 
that at least at the initial stages of the tube formation its temperature decreases exponen- 
tially, i.e. very sharply. The same conclusion can be made for the pressure and the matter 
density in the tube. Let us show it. 
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Assuming that the relation (B.15) holds not only for r = 0, but also for small r close to 
zero, we obtain 

Pintit) = PintiO) exp(-7t), (B.18) 

i.e. the inner pressure also decreases exponentially, but with the different exponential factor. 
Further, assuming that the heat capacity of a molecule in the tube is 

c=l- (B.19) 

where i is a number of molecule's degrees of freedom, and substituting the expressions (B.17) 
and (B.18) into the tube's equation of state, we derive the expression for the matter density in 
the tube 



,,. At* PintiO) 



1 - - ) 7t 



(B.20) 



Taking into account that z ^ 3, we can see that the density decreases exponentially just like 
the temperature (B.17) and pressure (B.18). 

B.5 Ideal photon channeling (without absorption) conditions 
inside the magnetic flux tubes 

Let us calculate the inner pressure pmt of the magnetic flux tube in the solar tachocline. In 
the framework of the standard model of the Sun the pressure in the tachocline zone is about 
~ 6- lO^'^ Pa, while the magnetic field strength reaches 400 T, according to our estimates (Fig. 9 
in the body text). Then from the hydrostatic condition by Parker (B.8) it follows that the inner 
pressure of the magnetic flux tube in the solar tachocline is equal to 

K„f~ 6- 10^2 -4- 10^° -6-10^2 [Pa], (B.21) 

which is comparable to the external pressure. 

However, the situation changes drastically within the framework of the "axion" model of the 
Sun. We pointed out earlier (A) that the values of thermodynamic parameters (temperature, 
pressure and plasma density) in the "axion" model are substantially smaller than the corre- 
sponding parameters in the standard model of the Sun. Let us suppose that the pressure in the 
tachocline zone in the "axion" model falls by an order of magnitude and is about ~ 10^^ Pa. In 
such case it is easy to see that the magnetic field as strong as ~ 500 T compensates the outer 
pressure almost completely, and the inner pressure 

Pi„i ~ 10^1 - O(10") ^ [Pa] (B.22) 

becomes ultralow. 

The result (B.22) means that the temperature and the matter density decrease together 
with the inner pressure, and the decrease is sharply exponential, because the exponential factor 
in (B.16) becomes very large 

7r = ^^ - 1 ^ oo. (B.23) 

Pint{T) 

Because of the fact that the density, pressure and temperature in the tube are ultralow, 
they virtually do not infiuence the radiation transport in these tubes at all. In other words, the 
Rosseland free paths for 7-quanta inside the tubes are so big that the Rosseland opacity and 
the absorption coefficients in (A. 32) tend to zero. The low refractivity (or high transparency) 
is achieved for the following limiting condition: 
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Pext ^ 1^. (B.24) 

The obtained results should not be considered as a proof, but rather as some trial estimates 
that validate the substantiation of an almost ideal (without absorption) photon channeling 
mechanism along the magnetic flux tubes. 

A major disadvantage of our reasoning is that, first of all, our estimates apply to the initial 
magnetic tube formation stages, and second, that it lacks the model calculations of the temper- 
ature, pressure and density in the framework of the "axion" model of the Sun, and consequently, 
there is no comparison with the corresponding parameters of the standard model. On the other 
hand, we believe that the idea of 7-quanta channeling along the magnetic flux tubes is so phys- 
ically natural and promising, that these disadvantages will be overcome in the near future in 
spite of the obvious serious difficulties. 
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